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ABSTRACT 


Few  underwater  acoustic  models  have  been  validated  in  the  Very  Low  Frequency  (VLF) 
regime.  The  purpose  here  is  to  explore  one  model’s  abilities  and  compare  it  widt  ocean  data  once 
it  is  ascertained  the  model  has  converged  to  its  final  solution.  A  comparative  analysis  of  data 
recorded  in  NATIVE  1  (the  Noise  and  Transmission  Loss  in  VLF  Environments)  Event  1  and  results 
generated  using  the  Seismo-Acoustic  Fast  field  Algoridun  for  Range-Independent  Environments 
(SAFARI)  is  die  focus  of  this  study,  widi  emidiasis  on  understanding  some  of  die  model’s  features. 
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I.  IMTRODDCTION 


In  submarine  warfare  there  is  a  need  for  accurate  acoustic 
models  both  for  the  prediction  of  how  well  a  given  sensor  will 
work  and  to  support  the  design  of  new  sensors.  These  acoustic 
models  need  to  be  verified,  in  the  acoustic  regime  of 
interest,  against  a  database  of  experimentally-measured  data 
representative  of  the  environment. 

The  verification  of  an  acoustic  model  requires  three 
primary  ingredients.  The  first  is  a  model  that  produces 
results  that  are  testable,  the  second  is  a  database  against 
which  to  compare  the  results  of  the  model,  and  third  is  a 
method  of  comparison;  whether  statistical  or  graphical.  One 
purpose  of  any  acoustic  model  is  to  simulate  the  environment 
with  sufficient  accuracy  so  results  can  be  compared  to 
measurements.  Limitations  on  the  model  input  parameters  have 
direct  ramifications  on  its  results.  Similarly  the  precision 
and  accuracy  of  the  equipment  used  in  recording  the  data  as 
well  as  any  perturbations  which  may  occur  during  the 
collection  of  the  data  must  be  taken  into  account  when 
determining  its  quality.  A  model  is  only  as  good  as  its 
theoretical  design  and  an  experiment  is  only  as  good  as  its 
design  and  equipment.  Once  both  results  are  compiled  they  can 
be  compared.  Then  a  determination  can  be  made  as  to  whether 
the  model  matches  with  the  features  of  the  measured  data. 
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The  objective  of  this  research  is  to  verify  the  Seismo 
Acoustic  Fast  field  Algorithm  for  Range- Independent 
environments  (SAFARI)  in  the  very-low-frequency  (VLF)  acoustic 
regime  (from  5  Hz  to  20  Hz)  by  comparison  with  the  data 
recorded  by  Marine  Physical  Laboratory  (MPL)  in  the  Noise  and 
Transmission  Loss  in  Very  Low  Frequency  Envi  ..:>nment  (NATIVE  l) 
experiment,  which  they  conducted  in  August  of  1990.  Chapter 
II  gives  a  thorough  background  on  the  NATIVE  1  experiment, 
including  a  description  of  the  site's  geoacoustic  and 
bathymetric  makeup  as  well  as  the  types  of  source  and  receiver 
used.  Chapter  III  goes  into  a  description  of  the  theory 
behind  SAFARI.  Chapter  IV  discusses  the  data  sources  used  in 
the  comparison.  Chapter  V  investigates  model  aliasing  and 
contour  offset  features.  Chapter  VI  presents  a  graphical 
comparison  of  model  results  and  data.  Chapter  VII  gives  the 
conclusions  reached  and  proposes  areas  for  further 
investigation . 
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IZ.  NATIVE  1 


A.  BACKGROUND 

In  August  of  1990,  the  NATIVE  1  experiments  were  conducted 
to  acquire  data  to  provide  insight  into  VLF-propagation 
effects.  Along  with  the  development  of  a  VLF  data  base  was 
the  effort  to  design  VLF  performance  prediction  models  which 
would  support  the  development  of  new  Antisubmarine  Warfare 
(ASW)  acoustic  sensors.  NATIVE  1  was  a  multi-laboratory, at- 
sea  exercise  conducted  on  Blakes  Plateau,  Figure  1 
[Ref.  1].  Two  sites  were  used,  one  in  deep  water  and 
one  in  shallow  water.  Five  events,  lasting  several  weeks  were 
conducted  using  a  variety  of  devices  and  equipment.  Data  from 
the  first  of  these  events  at  the  deep  water  site  will  be  used. 
Figure  2  [Ref.  2]  shows  the  locations  of  the  initial 
entry  and  recovery  points  of  the  receivers.  Typical  Swallow 
Float  deployments  consist  of  forming  an  equilateral  triangle 
pattern  with  three  floats,  at  known  positions,  spaced  about 
6.5  km  between  floats.  This  provides  a  reference  for  float 
localization.  A  15-lb  anchor  on  a  3-m  tether  replaces  the 
ballast.  The  remaining  floats  are  deployed  in  the  center  of 
the  triangle.  Figure  3  shows  the  path  of  the  source  over  the 
day  and  a  half  comprising  the  recording  time.  The  event 
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Figure  l  Native  l  site  "D**  in  3400  m  of  water  and  "S”  in 
1000  m  of  water. [Ref.  l:p.  2] 

starts  with  the  source  well  north,  approximately  12  km,  of  the 
center  of  the  sensor  site.  It  maintains  a  track  of 
approximately  193  degrees  past  the  receivers  and  continues 
south  for  another  50  km,  and  then  turns  around  on  the 
reciprocal  heading  finishing  approximately  6  km  north  of  the 
site. 

B.  IZTRBMBLT  LOW  FREQUENCY  SOURCE 

In  this  event  a  calibrated  extremely  low  frequency  (ELF) 
source  projector  was  towed  by  the  research  vessel  (R/V)  Gyre. 
A  precise  positional  log  (Appendix  A)  was  maintained  from 


4 


Flgur*  2  Plan  view  of  entry  and  recovery  positions  of 
Swallow  floats.  Event  1,  NATIVE  1,  (34 *N,  74 ‘W)  [Ref  2:p.  2] 


Miiwta 

,  Figure  3  Research  vessel  Gyre  track  during  Event  l,  NATIVE 

1  (31*N,  74*W)  [Ref.  2;p.  2] 

•  multiple  navigational  sources.  These  included  Global 

Positioning  Systems  (GPS) ,  LORAN,  and  Transit  Satellite 
Navigation  (Transit  Sat  Nav) .  Towing  speed  was  4.0  ±  0.5 
knots  and  the  source  depth  was  122  ±  7  meters.  The  source 
projected  two  continuous  wave  frequencies  simultaneously;  7  Hz 
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and  10  Hz  from  1630  to  2230  on  4  August  and  0030  to  0400  on  5 
August,  and  10  Hz  and  16  Hz  from  0800  to  1130  on  5  August.  At 
7  Hz,  the  average  source  level  was  164  ±  2  dB//jiiPa,  at  10  Hz 
the  average  source  level  was  166  ±  2  dB/^Pa,  and  at  16  Hz  the 
average  source  level  was  161  ±  2  dB//xPa  [Ref.  3]. 
(All  estimates  were  made  from  visual  Inspection  of  the  data.) 
The  positional  accuracy  of  the  ship  was  source  dependent  with 
the  GPS  system  having  the  highest  degree  of  accuracy,  about  ± 
60  meters.  Transit  Sat  Nav  was  next  with  an  accuracy  of  ±  100 
meters  and  LORAN  was  accurate  to  within  ±0.5  nautical  miles. 
The  times  each  of  these  systems  were  used  on  the  above  dates 
are  highlighted  In  Appendix  A,  along  with  the  actual  track 
data. 


"The  Transit  Sat  Nav  Is  dependent  on  satellite 
elevation  above  the  horizon,  time  from  satellite  fix 
to  position  solution  and  dead  reckoning  parameters 
(course  and  speed)  used  for  updates.  Therefore  I 
estimate  Transit  Sat  Nav  Is  only  accurate  In  general, 
to  within  O.S-l.O  nautical  miles.  LORAN  appears  to  be 
accurate  to  within  0.5  nautical  miles  with  better 
repeatability.  For  this  reason  there  Is  a  constant 
offset  or  position  error  between  GPS  and  LORAN  of  <• 

0.021  degrees  longitude  which  should  be  added  to  LORAN 

longitudes.  This  would  move  all  LORAN  positions 

further  to  the  west.  This  correction  has  not  been 

applied  to  LORAN  positions  In  enclosure (2)  (Appendix 

A)  which  Is  taken  from  the  NADC  navigation  « 

rectification  products.  It  also  will  not  affect 

ranges  along  the  track  from  source  to  receiver  due  to 

the  north-south  orientation  of  the  Event  l  track."  * 

[Ref.  4] 
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C.  BllVIROIlMEliTAL  DATA 

The  exercise  area  is  located  off  Jacksonville  Florida, 
east  of  the  Gulf  Stream  in  the  waters  commonly  called  the 
Sargasso  Sea.  At  the  end  of  Event  1,  41  bathythermographs  and 
one  conductivity,  depth,  temperature  (CDT)  cast  were  made 
(Figure  4) .  Sample  values  of  those  results  [Ref.  l:pp.  Bl- 
B23]  are  in  Appendix  B,  and  are  typical  of  the  water  masses 
which  commonly  make  up  this  region.  The  temperature 


meoiATUc  OO  (f) 

Figure  4  Site  "D"  temperature  profiles  for  (a)AXBTs,  (b)XBTs 
[Ref.  l:p.  22] 


variability  was  minimal  with  one  exception  located  in 
Gulf  stream  Haters  well  to  the  west  of  the  track.  Province 
contours  (Figure  5)  were  generated  and  representative  values 
of  each  province  are  listed  in  Appendix  B.  Those  profiles 
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Kuiw:  I  ;u4M»  nri»  riioKiu;  piiM>ini-i::i 


-7*  -75  -V4 

UMCI1UK  1DI1C> 


Flgur*  S  site  "O'*  sound  speed  provinces.  [Ref .  l:p.  23] 

that  are  generally  cooler  have  lower  province  numbers.  Table 
I  shows  the  province  profile  number  and  range  from  site  "D" 
along  the  track  of  R/V  Gyre  during  Event  l. 

Seasonal  heating  in  summertime  dominates  the  surface  water 
while  the  seasonal  thermocline  extends  to  about 
100  m.  From  100  to  500  m  is  the  regime  called  *18  degree' 
water  mass  where  the  temperature  gradient  is  small  producing 
an  isovelocity  water.  Below  500  m  to  any  where  between  2500 

to  3000  m  North  Atlantic  Central  Water  with  a  deep  sound 

Table  Z  RANGE  FROM  SITE  "D**  AT  WHICH  THE  SOURCE  ENTERED  A 
PROVINCE  PROFILE  NUMBER  ALONG  EVENT  1  TRACK. [Ref  lip. 24] 
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channel  axis  at  about  1200  meters.  Below  this  is  the  very 
stable  North  Atlantic  Deep  and  Bottom  water  mass. 

D.  6E0AC0U8TIC  BOTTOM 

The  bathymetric  track  (Figure  6)  was  recorded  by  R/V  Gyre 


40} 

during  Event  1.  The  bottom  varied  from  a  minimum  of  3159  m 
near  the  VEDABS  1  location  (center  of  the  site)  to  a  maximum 
at  3688  m.  This  yields  an  approximate  slope  from  0.002  to 
0.001  for  the  bottom  and,  as  such,  is  considered  horizontal 
and  flat.  In  Figure  7  representative  samples  of  the 
geoacoustics  of  the  bottom  sediments  from  [Ref.  5] 
were  generated  along  the  R/V  Gyre  track  using  the  model  found 
in  the  front  of  Appendix  C. [Ref  l:p.  45] 

E.  SHALLOW  FLOATS 

The  following  discussion  is  adapted  from  D'Spain  et.al. 
[Ref.  6].  A  general  description  of  the  Swallow  Float 
is  followed  by  a  discourse  on  the  data  acquisition  process, 
including  a  discussion  of  the  hydrophone  and  the  electronic 
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BATHYMETRY 

DEPTH  TO  ACOUSTIC  BASEMENT 


RANGE (KM) 


Figur*  7  site  "D"  Bearing  000°T  Water  Depth  and  Depth  to 
Acoustic  Basement  [Ref.  5:p.  B-‘24] 


circuit,  and  finally  a  discussion  on  the  hydrophones 
autospectra . 

1.  Description 

Swallow  Floats  are  independent,  freely-drifting, 
neutrally-buoyant  sensors  which  can  be  ballasted  to  float  at 
any  desired  depth  (±  100  m) .  Figure  8  shows  a  schematic  of  a 
typical  Swallow  Float.  Mounted  within  the  0.432m  diameter 
Benthos  glass  shell  are  three  orthogonally  mounted  geophones, 
a  compass,  the  electronics  with  power  supply  (not  shown) ,  and 
cassette  tape  recorder  (not  shown) .  A  polyethylene  hard  hat 
covers  the  glass  shell's  exterior  with  an  Ocean  Applied 
Research  (OAR)  xenon  flasher  attached  to  one  side  and  an  OAR 
radio  beacon  on  the  other.  These  are  used  for  detection  and 
localization  during  recovery. 


To  prevent  distortion  of  the  compass  readings,  the 
compass  is  mounted  on  the  bottom  of  the  shell  at  its  maximum 
‘iistance  from  the  magnetic  fields  of  the  geophones,  mounted  on 
the  top  of  the  shell.  The  power  supply  provides  a  maximum  of 
200  hours  of  data  acquisition  and  25  days  of  acoustic- 
transponder  life.  Duration  of  Swallow  float  deployments  is 
not  limited  by  power  supply  requirements  but  by  the  capacity 
the  17-megabyte  cassette  tape  recorder. 
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An  infrasonic  hydrophone  is  mounted  outside  the  glass 
shell  on  the  base  plate  just  below  the  compass,  approximately 
0.5  m  below  the  geophones.  Suspended  below  the  base  plate  is 
the  expendable  ballast.  The  ballast  can  be  acoustically 
commanded  to  release  and  has  two  independent  backup  timed 
release  circuits.  The  ballast  weight  is  accurately  determined 
prior  to  float  deployment  to  achieve  the  desired  depth  of 
deployment.  A  change  in  ballast  of  1  gm  results  in  depth 
changes  of  7  to  10  m.  The  ballast  fits  snugly  against  the 
base  plate  in  order  to  reduce  its  rocking  motion  relative  to 
the  float;  herein  Figure  8  is  distorted. 

Each  float  is  also  equipped  with  an  International 
Transducer  Corporation  (ITC)  localization  hydrophone  which  is 
suspended  1.83  m  below  the  glass  shell.  This  hydrophone  has 
a  source  level  of  191  dB//lMPa  at  1  m  and  a  receiving 
sensitivity  of  -180  dB//lV/MPa  at  8  kHz.  Under  microprocessor 
control  and  synchronized  prior  to  deployment,  the  floats  send 
8-kHz  localization  pulses  of  lO-ms  duration.  This  allows  for 
post  deployment  determination  of  depth  and  location  during 
deployment.  This  hydrophone  and  the  associated  circuitry  also 
allow  for  two  types  of  acoustic  commands,  uniquely  coded  for 
each  float,  to  be  issued.  The  first  is  a  50-ms  8 -kHz  pulse 
used  for  ship  to  float  slant  range  calculation  and  the  second 
activates  the  ballast  release  circuitry  and  issues  a  500-ms  8- 
kHz  pulse  and  is  used  at  the  end  of  the  experiment  for  float 
recall. 
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2.  Dftta  Aoquisitlon 


r 


% 


The  Swallow  Float  data  acquisition  system  (Figure  9) 
records  the  acoustic  pressure  and  particle  motion  at  the 
sensors.  This  discussion  is  restricted  to  the  acoustic 
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Figure  9  swallow  Float  data  acquisition  system  [Ref.  6:p. 
31] 


pressure  only.  An  OAS  model  E-4SD  hydrophone  is  the  pressure 
sensor  which  converts  acoustic  and  nonacoustic  pressure 
fluctuations  into  voltage  signals  with  a  sensitivity  of  -182 
dB//  1  V//iiPa.  The  hydrophone  has  a  frequency  response  which 
is  flat  (±  1  dB)  from  0  to  5kHz  [Ref  6:p.  8].  The  signal  then 
passes  in  order  through: 
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•  two  prewhitening  RC  filters 

•  fixed  80  dB  gain  amplifier 

•  three  highpass  RC  filters  and  four  lowpass  RC  filters 
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•  automatic  gain  control  (AGC)  amplifier 

•  anti*-aliasing  filter 

•  A/D  converter 

•  data  buffer  with  44  seconds  of  data 

•  cassette  tape  recorder  with  1  second  of  write  time  for 
each  acoustic  record 

Each  acoustic  record  is  45  s  long  (44  seconds  of  recorded 
data  and  1  second  of  write  time) . 

3 .  Pressure  Autospeetra 

From  the  recorded  data,  the  RMS  pressure  at  each  float 
was  calculated.  These  RMS  time  series  were  corrected  for  AGC 
gain  but  were  otherwise  uncalibrated.  The  vertical  axis  of 
the  recorded  data  is  in  units  of  RMS  volts  at  the  A/D 
converter  (Figure  9) . 

"Each  RMS  level  value  results  from  taking  the  square  root 
of  the  average  of  the  squared  amplitude  levels  over  a 
period  of  five  seconds  (or  250  points  since  the  data 
sampling  rate  is  50  Hz) .  Since  each  record  is  2250  points 
(45  seconds)  long,  then  nine  RMS  power  values  are 
calculated  for  each  record.  The  RMS  value  for  the  last 
five-second  period  of  each  record  includes  one  second  of 
zeros.  The  one  second  of  zeros  represents  the  time  during 
which  no  data  is  sampled  while  data  in  a  temporary  buffer 
is  being  written  to  cassette  tape."  [Ref.  2:p.  21] 


"Spectral  estimates  with  two  different  frequency 
resolutions  are  calculated.  Those  with  a  98  mHz 
resolution  are  made  by  dividing  40.96  seconds  in  a  data 
record  (obtained  after  skipping  the  first  3  seconds  of 
data)  into  seven  10.24-second-long  segments,  with  a  50 
percent  overlap  between  segments.  The  segments  are  fast 
Fourier  transformed  after  using  a  Kaiser-Bessel  window  of 
a*2.5  to  window  the  data.  The  resulting  seven  spectra 
within  a  record  are  incoherently  averaged  together.  The 
averaged  spectra  for  four  consecutive  records  (equal  to 
three  minutes  of  data)  are  then  incoherently  averaged  in 
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order  to  further  reduce  the  variance  of  the  spectral 
estinates.  For  spectral  estimates  with  a  24  mHz 
resolution  ,  the  segment  length  which  is  transformed  is 
quadrupled  to  40.96  seconds  and  the  spectra  for  twenty 
consecutive  records  are  incoherently  averaged. " [Ref .  2:p. 
26] 

The  autospectra  are  calculated  in  dB  //  l(/xPa)2/Hz.  Records 
used  in  the  analysis  were  all  three  minute  averaged  data, 
using  four  45  second  records  for  each  data  point. 
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III.  MXTBEMXTICAL  MODEL 


A.  BACXORODllD 

In  the  attempt  to  model  the  complex  environment  found  in 
the  ocean,  there  have  been  a  variety  of  techniques  used  to 
determine  the  solution  to  the  wave  equation,  each  involving 
its  own  simplifying  assumptions  and  approximations.  One 
subset  of  these  techniques  is  the  method  using  separation  of 
variables.  Two  approaches  fall  into  this  general  category: 
the  normal  mode  and  the  fast  field  approach.  It  is  the  latter 
that  is  Implemented  in  SAFARI  [Ref.  7].  The  following 
sections  provide  a  summary  of  the  wavefield  technique  found  in 
Schmidt  [Ref  7].  In  a  horizontally  stratified  environment, 
including  fluid  and  solid  layers,  this  technique  yields,  an 
exact  solution  to  the  linearized  wave  equation, 

(V* — i — --^)T(r,  z,  t)  »F,(r,z,  t)  3.1 

c(r,z)>at* 

The  following  sections  develop  the  basic  mathematical 
approach,  discuss  environmental  and  analytical  models, 
boundary  conditions,  and  numerical  techniques  including,  the 
global  matrix  approach,  wavenumber  integration,  aliasing,  and 
contour  offsetting. 

B.  MATHBIATXCAL  MODEL 

As  early  as  1948  [Ref.  8],  the  subject  of  wave 
propagation  in  a  stratified  medium  was  being  studied  in  planar 
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two-  and  three-layer  environmental  models.  The  technique 
applies  a  series  of  integral  transforms  to  the  linearized 
Helmholtz  wave  equation,  Eq.  3.1,  reducing  a  partial 
differential  equation  in  time  and  three  spatial  dimensions 
into  a  series  of  ordinary  differential  equations  in  the  depth 
coordinate.  These  ordinary  differential  equations  can  be 
solved  analytically  in  terms  of  unknown  amplitudes  by  matching 
boundary  conditions  at  the  interfaces.  For  a  cylindrical 
coordinate  system  {r,e,z},  the  model  solves  for  an  axial 
symmetric,  or  planar,  wavefield  independent  of  6.  The 
isotropic  media  can  be  expressed  in  terms  of  scalar  potentials 
i'(r,z,t)  which  satisfy  Equation  3.1,  where  c(r,z)  is  the  wave 
speed  and  F,(r,z,t)  is  a  forcing  term  representing  the  source. 
The  first  of  two  integral  transforms  is  next  applied.  The 
forward  Fourier  transform 


(a) 

2it  J 

F(t)  =]■/(«)  ib) 


3.2 


gives  the  frequency  domain  wave  equation 

(^*kjl(r,z))^{r,z,ia)  3.3 

where  u  is  the  angular  frequency  and  k.(r,z}  is  the  medium 
wave  number, 

Jc,(r,r)  »«/c(r,z)  3.4 
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For  stratified  range  Independent  environments  with  a  point 
source  present.  Equation  3.3  can  be  simplified  to 


2iir 


In  the  following,  the  u  dependence  will  be  suppressed. 
HanXel  transform  pair  are  utilized 


3.5 

The 


gik)  =jG{z)  Jo(kz)  zdr  (a) 
0 

G(z)=fg(k)Jg(kz)kdk  {b) 


where  k  is  the  horizontal  wavenumber,  and  the  forward 
transform,  Eq  3.6(a) ,  is  applied  to  obtain  the  depth-separated 
wave  equation  which  is  an  ordinary  differential  equation  in 
depth  only.  The  solution  to  this  equation  is  the  sum  of  a 
particular  solution to  Eq.  3.5,  and  any  linear  combination 
of  the  two  independent  solutions  to  the  homogeneous  equation. 
It  is  commonly  called  the  depth-dependent  Green's  function 
whose  coefficients  are  arbitrary  and  are  determined  by  the 
boundary  conditions.  Once  the  coefficients  are  found,  the 
inverse  Hankel  transform  [Eq.  3.6(b)]  yields  a  wavefield  at 
any  range  r  at  the  angular  frequency  w.  The  inverse  Fourier 
transform  yields  the  time  response  in  the  same  fashion. 
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C.  BOUMDARY  CONDITIONS 


4 


1 .  Environmental  Model 

In  SAFARI,  the  environment  is  represented  as  a 
horizontally  stratified  medium.  Figure  10,  with  all  interfaces 
planar  and  parallel.  The  properties  of  the  layers  are  range 
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Tigure  10  Horizontally  stratified  environment  [Ref.  7:p.  5] 


independent.  Solid  layers  are  assumed  to  be  homogeneous  and 
isotropic  with  Lam6  constants  X  and  ju  and  density  p.  The 
compressional  wave  speed  is 

Cg=^/ln^iiT7p  3  •  7 
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and  the  shear  speed  is 


3.8 


Physically  realistic  materials  must  satisfy  the  condition 


c.<v^r7?c^ 


3.9 
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so  that  the  bulk  modulus,  K=  X+2/3/x,  be  positive,  and  the 
solid  medium  have  positive  compressibility.  The  density  in 
fluid  layers  must  be  constant  while  the  nonvanishing  Lamd 
constant  X  is  allowed  to  vary  with  depth  according  to: 

V”  ■  =az*b  3.10 

X(2) 

corresponding  to  the  following  depth  dependence  for  the  sound 
speed,  Eq  3.7, 

c^iz)  =V1/  (P  (az*Jb) )  3.11 


2 .  Analytical  Models 

The  analytical  field  representation  for  the  media 
included  in  SAFARI  are  for  a  homogeneous  fluid,  inhomogeneous 
fluid,  and  solid  representations,  and  are  governed  by  the 
equation  of  motion  for  each  media.  The  focus  of  this  section 
is  not  on  the  details  of  each  derivation  but  rather  on  the 
fact  that  the  wavefield  at  each  Interface  has  two  integral 
representations  one  each  from  the  upper  and  lower  layers. 

For  the  homogenous  fluid,  the  solution  to  the  depth- 
dependent  Greens  function  for  the  potential  is  of  the  form 


♦  (2,2)*!  [A-e-^’*A*e*’J^{kr)  kdk  3.12 

0 
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where  k„(2)»h„  (a  constant),  a  (k)  -yjk^-hl  ^  and  no  source  is 

present.  A'  and  A^  are  the  unknown  amplitudes  of  the  downgoing 
and  upgoing  waves.  If  an  omnidirectional  point  source  is 
present  the  particular  solution  is  added  to  the  homogeneous 
solution  and  is  of  the  form 

*  (r,  r)  f  ^  Jp (Jcr) kdk  3.13 

4tC<  ft 

where  is  the  source  strength  and  z,  is  the  source  depth.  The 
associated  displacement  and  stress  are 

w{z,  z)  fsign(z-Zg)  e'***’*''^^  ikr)  kdk  3 . 14 

4tc  J 


Oggiz,  z) 


i  a 

0 


(itr)  kdk 


3.15 


where 


I 


aigniz-Zg) 


l,Z-Zg>0, 
0,  z-z,=0, 
-1,  Z-Z,<0 


Since  the  ocean  is  not  a  homogeneous  media  and  the 


sound  velocity  varies  with  depth,  a  single  layer  ocean  wave 
guide  doesn't  properly  reflect  acoustic  propagation.  However, 
if  a  range  independent  wave  guide  is  represented  by  an 
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number  of  homogeneous  layers,  a  numerical  solution  based  on 
their  field  representation  will  converge  toward  the  correct 
solution.  It  can  be  shown  when  modelling  a  deep  water  site 
there  must  be  a  velocity  profile  incorporated,  as  is  the  case 
here,  moreover,  to  ensure  a  satisfactory  convergence  the 
isovelocity  layers  must  be  less  than  one  quarter  wavelength 
thick.  Once  again  it  can  be  shown  the  solution  to  the  depth 
dependent  Greens  function,  based  on  c(z),  from  Eg.  3.11,  and 
not  an  isovelocity  constant,  can  be  obtained  by  using  a 
variable  transformation  in  which 


When  substituted  into  Eg.  3.5  the  following  is  obtained: 


(^-C)*(C)*0 

av 


3.17 


This  yields  the  integral  form  of  the  potential  with  a  source, 
®  (r,  z)  =--^  fjo 

_ (Ai  (C,)  ♦iBi  (C,) )  Ai  (C) _ 

Ai'(C,)  Ui(C,)  ♦iBi(C.))-Ai(C,)  (Ai'(C,)  ♦iBi'(C,)) 

3  •  18 

X,  ,a(z-z,)iO, 

2c-i/5Ai(C.)  Ui(C)+iBi(C)) 

Ai'tC.)  <Ai{C.)  +ifli(C.))-Ai(C.)  (Ai'(C.)  +iBi'(C.)) 

,a(z-z,)  iO 
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By  expressing  the  displacements  due  to  both  normal 
stress  and  tangential  stress  forces  in  terms  of  scalar 
potentials  the  solution  to  the  depth  dependent  Greens  function 

for  a  solid  are  of  the^form, 

«(r,r)=| tA-e-**+A*e**]  J-o(;cr);cdic  3.19 

0 

'9{r,z)  •f  iB-e'^‘*B*e^*]Jo{kz)dk  3.20 


The  normal  stress  and  tangential  stress  are  solved  for 
in  the  same  fashion  as  in  the  homogeneous  case.  Here  0  is  the 
wavenumber  for  tangential  stress  analogous  to  the  normal 
stress  wavenumber  a  in  Eq.  3.12. 

3 .  Msteliiag  Boundaries 

The  boundary  conditions  are  matched  at  each  interface 
dependent  upon  the  type  of  interface  involved.  The  following 
list  summarizes  the  conditions: 

•  For  an  interface  separating  two  fluid  layers,  the  vertical 

displacement  w  and  the  normal  stress  have  to  be 

continuous.  If  one  of  the  media  is  a  vacuxim,  the  normal 
stress  must  vanish. 

•  For  an  interface  separating  a  solid  and  a  vacuum,  both 
and  must  vanish,  whereas  a  solid-fluid  interface 
requires  that  w  and  must  be  continuous  and  a„  must 
vanish. 

•  For  an  interface  separating  two  solid  media,  w,u,a„,  and 
a„  must  all  be  continuous. 

"Since  the  boundary  conditions  have  to  be  satisfied  at  all 
ranges  r,  it  is  obvious  that  they  must  be  satisfied  by  the 
kernels  in  the  Integral  representations  as  well.  The 


23 


radiation  conditions  for  z-*±»  together  with  the  conditions 
to  be  satisfied  at  all  interfaces  simultaneously,  lead  to 
a  linear  system  of  equations  in  the  unknown  kernel 
coefficients  and  B*.  In  principle,  this  system 

has  to  be  solved  for  all  values  of  the  horizontal 
wavenumber  k,  and  the  total  field  can  then  be  determined 
by  evaluating  the  inverse  transforms.  Except  for  a  few 
trivial  cases,  however,  both  the  solution  of  the  linear 
system  of  equations  and  the  evaluation  of  the  inverse 
transforms  have  to  be  done  numerically,  requiring 
truncation  and  discretization  of  the  horizontal  wavenumber 
axis."  [Ref.  7:pp.  17-18] 

D.  VDMSRZCAL  TBCBMIQUB8 

The  global  matrix  approach  and  wavenumber  integration  are 
now  discussed. 

1.  Global  Matrix  Approach 

The  global  matrix  approach  is  a  more  direct  and 

computationally  more  efficient  technique  in  assembling  the 

local  solution  to  the  depth  dependent  Greens  function,  for 

each  layer,  in  terms  of  its  source  contribution  and  unknown 

scalar  potentials,  with  all  other  local  solutions.  The 

boundary  conditions  determine  a  system  of  equations  in  the 

Hankel  transforms  of  the  potentials  to  be  satisfied  at  each 

interface.  In  similar  fashion  to  finite  element  model 

programs  [Ref.  9]  these  local  system  of  equations  are 

mapped  into  a  global  set  of  equations  defining  the  boundary 

conditions  at  all  interfaces  simultaneously. 

"The  resulting  global  coefficient  matrix  is  organized  to 
be  block-bidlagonal  and  diagonally  dominant  in  close 
analogy  to  the  global  stiffness  matrix  arising  from  the 
finite-  element  method.  Unconditionally  nvimerically- 
stable  solutions  are  determined  efficiently  by  gaussian 
elimination,  yielding  the  field  in  all  layers 
simultaneously."  [Ref  7:p.  29] 
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The  nechanics  of  the  technique's  derivation  of  the  solution  is 
secondary  to  the  understanding  that  unconditionally  stable 
solutions  are  obtained,  and  that  the  technique  is  stable  for 
any  machine  precision.  Herein,  SAFARI  is  an  excellent  choice 
for  depth  range  contouring  of  single  frequency  transmission 
loss,  the  focus  of  this  study. 

2 .  Wavenumber  Integration 

Numerical  techniques  which  evaluate  the  depth 
dependent  Green's  function  by  the  inverse  Hankel  transform, 

m 

G(T,2)*jg(k,z)J^{kr)kdk  3.6(b) 

0 

derive  truncated  solutions  to  the  infinite  integration 
interval.  Here  the  subscript  m  equals  0  except  for  the  shear 
stress  and  horizontal  displacement  where  m  equals  l.  Because 
the  interval  on  the  wavenumber  k  goes  from  0  to  »  the  solution 
can  be  done  to  any  degree  of  accuracy  within  the  computing 
limits  of  the  machine.  In  a  lossless  media,  waveguide 
propagation  produces  normal  modes  and  interface  waves  which 
show  up  in  the  wavenumber  integration  as  poles.  As  the 
Integrand  varies  over  k  these  plot  in  the  complex  plane  on  the 
real  axis.  For  lossy  media,  these  poles  move  out  into  the 
complex  plane,  remaining  close  to  the  real  axis.  In  cases 
where  solid  layers  underlie  fluid  layers  care  must  be  taken  to 
ensure  the  maximum  medium  wavenumber  truncation  point  is 
selected  beyond  any  significant  contributions  from  waves 
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propagating  in  those  solid  layers.  Here,  evanescent  surface, 
interface  or  plate  inodes  may  exist  causing  poles  to  appear  in 
the  integration  plots.  Typically  20%  above  a  maximum  medium 
wavenumber  is  a  good  value  for  truncation  of  the  wavenumber 
range . 

Once  the  truncation  point  is  selected  two  factors 
complicate  the  remaining  solution.  The  first  is  the  existence 
of  normal  modes  and  interface  waves  which  generate  poles  in 
the  kernel  of  the  transform.  The  second  is  the  Bessel 
function  J„(kr)  which  may  vary  rapidly  in  k,  the  horizontal 
wavenumber,  for  long  ranges  in  r.  A  very  fine  sampling 
interval  has  to  be  chosen  to  avoid  aliasing,  addressed  below. 

DiNapoli  and  Deavenport  [Ref.  10}  show,  for 
the  ranges  of  this  study,  the  Bessel  function, 

)  3.21 

can  be  represented  by  the  outgoing  Hankel  function 


Urn  (kz) 


2  (*»-!>  ■I"’ 


3.22 


(The  incoming  Hankel  function  is  needed  only  for  evaluation  in 
the  near  field  within  a  few  wavelengths  of  the  source.)  This 
yields  an  Inverse  Hankel  transform  integral 


Hm*h-kn 

I  • 


z)  ’/ke'^^dk. 


3.23 


« 


r 
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The  waveniunber  space  having  been  truncated  is  divided 
equidistantly  in  SAFARI  as  follows: 

iCj=ic;„in+lAic,  2=0,1.  ..  (W-1)  3.24 

where  M  is  the  total  number  of  sampling  points.  The  range 

axis  is  divided  accordingly  as 

j=0.1...(M-l)  3.25 

where  the  range  step  Ar  is  governed  by  the  relationship 

^TAk-2n/M  3.26 


and  M  is  an  integral  power  of  2.  The  discrete  approximation 
of  Eg  3.23  is  finally  obtained. 


Pi 

3.27 


The  sximmation  can  be  performed  by  means  of  a  fast  Fourier 
transform  (FFT)  generating  the  field  at  all  M  ranges,  in  Eg. 
3.25,  simultaneously. 

3.  Aliasing  and  Contour  offset 

Regarding  aliasing,  the  undersampling  in  one  domain 
causing  wrap  around  in  the  other  domain  of  a  transform.  Eg. 
3.27  does  not  yield  G(r,z)  but  a  summation  from  the  lower 
limit  to  the  upper  limit  of  the  interval.  The  range  from  zero 
to  the  lower  limit  and  from  the  upper  limit  to  infinity  will 
cause  some  aliasing  in  the  wavenumber.  Over  the  interval 
[r^,  r^  4-  R] ,  where  the  range  window  width  R-MAr,  the  maximum 
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useable  range  is  always  R»MAr,  independent  of  the  choice  of 


For  ranges  larger  than  r^+R  aliasing  can  be  reduced 
by  choosing  the  window  R  so  large  that  the  signal  is  known  to 
dissipate  within  the  window.  For  a  fixed  wave  number  interval 
this  would  require  a  large  number  of  wavenumber  samples. 

According  to  Cauchy's  theorem  the  integrals  in  the 
complex  plane  between  two  points  are  invariant  to  a  change  in 
the  integration  contour,  that  is  by  moving  the  integration 
contour  away  from  the  pole,  offsetting  it,  the  integrands  can 
be  smoothed  and  aliasing  avoided.  In  SAFARI  this  is  done  by 
replacing  k  with  ^»k-(-ie  where  e  is  a  vertical  offset  whose 
contribution  to  the  total  integral  is  insignificant  when 
compared  to  that  of  k  ie.  6«k^-k^.  Eq.  3.27  is  replaced  by 


G(z,z) 


N  2xr* 


2  2 


* jg(k,  3.28 


where  C  is  the  contour  shown  in  Figure  11.  For  ranges  larger 
than  r^  +  R,  a  factor  of  e"**  is  introduced  to  the  FFT  and  the 
signals  will  be  attenuated,  similarly  for  ranges  smaller  than 
signals  wrapped  around  will  be  amplified  by  e'*^.  No  change 
occurs  in  the  maximum  range  r_„»R,  from  the  real  axis 
integration.  Corresponding  to  an  attenuation  of  60  dB  of  the 
wrap  around  signal  the  value  of  e  is  formulated  as 

**  Rloge"  2x  (W-1)  loge 
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rigurtt  ll  complex  Integration  Contour  [Ref.  7:p.  33] 

which  even  for  a  small  nvimber  of  sampling  points  ensures  e  « 
Paramount  to  this  numerical  Integration  scheme  is 
that  the  kernel  must  vary  smoothly  between  sampling  points. 
If  the  Integration  is  truncated  at  a  point  where  the  kernels 
amplitude  is  significant,  the  discontinuity  will  give  rise  to 
aliasing.  To  avoid  this,  a  first  order  Hermite  polynomial  is 
employed  and  the  kernel  and  its  first  derivative  become 
continuous  at  the  truncation  point. 
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IV.  DATA  SOURCES  AND  ACOUSTIC  MODELING  INPUTS 
The  recorded  transmission  loss  data  sent  by  Marine 
Physical  Laboratory  (MPL)  of  Swallow  Float  0  is  described 
along  with  the  input  files  used  in  running  the  SAFARI  model. 

A.  SWALLOW  FLOAT  DATA 

The  data  [Ref.  11]  received  from  MPL  were  12 

files  of  Swallow  Float  0  data  recorded  during  Event  l  of  the 

NATIVE  1  experiment.  It  was  received  via  electronic  mail  and 

transferred  to  5.25**  high  density  floppy  disks.  The  files 

contained  two  columns  of  information.  The  first  column 

contained  the  R/V  Gyre-to-buoy  range  in  kilometers  and  the 

second  column  contained  either  the  percentage  of  clipped  data 

points  (four  files)  or  the  transmission  loss  (eight  files)  in 

units  of  dB  referenced  to  a  source  with  pressure  amplitude  of 

1  mPs  St  1  meter.  Reference  11  received  from  MPL  goes  into 

the  scheme  for  defining  the  data:  TLO. [1,2, 3,4] . [7,10,16] 

references  the  Transmission  Loss  Float  0,  the  region  number 

and  the  frequency  for  the  region.  The  regions  are. 

Region  l:  The  Gyre  approached  the  float  from  the  north- 
northeast  near  the  beginning  of  Event  1  broadcasting  7  and 
10  Hz. 

Region  2:  The  Gyre  passed  over  the  float  heading  south- 
southwest  broadcasting  7  and  10  Hz. 

Region  3:  After  the  SUS  charge  events  as  the  Gyre  continued 
to  the  south-southwest  broadcasting  7  and  10  Hz. 
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Region  4:  The  Gyre  approached  the  float  from  the  south- 
southwest  broadcasting  10  and  16  Hz. 

With  regard  to  the  navigational  data  and  the  Gyre-to-f  loat 
range,  the  5  minute  log  in  Appendix  A  was  interpolated  every 
45  seconds  to  get  a  position  for  each  record  of  data. 

As  of  the  writing  of  this  paper,  there  is  a  question  as  to 
whether  or  not  a  range  offset  of  approximately  0.5  km  occurs 
in  the  navigational  data. 

Unless  otherwise  annotated,  all  Swallow  Float  data 
presented  in  the  analysis  section  has  less  than  5  percent  of 
clipped  time  series  points  in  the  hydrophone  channel.  Where 
the  percentage  of  clipped  points  exceeded  5  percent  the  record 
was  deleted.  This  was  done  using  a  DOS  text  editor. 

B.  BJUPARZ  MODEL  AMD  INPUTS 

The  SAFARI  model  was  installed  at  the  Fleet  Numerical 
Oceanography  Center  in  Monterey,  California  in  a  UNIX 
environment  on  an  HP  workstation.  It  was  supported  by  DISSPLA 
graphics,  and  executed  the  Fortran  code  on  an  F77  Fortran 
compiler.  Upon  satisfactory  duplication  of  the  results  for  the 
test  files  in  Ref.  7,  and  an  independent  two-layer  split  fluid 
case  verified  by  alternate  means,  analysis  proceeded  with  the 
model. 

The  water  column  data  from  province  4  in  Appendix  B  and 
the  geoacoustic  model  data  from  Appendix  C  was  used  to  produce 
a  43-layer  water-sediment  input  file  in  the  standard  SAFARI 
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FIP  format.  The  FIP  module  generates  single-frequency 
calculation  of  the  total  wavefield  for  a  single  point  source 
whose  source  strength  is  normalized  to  yield  a  pressure  of 
1  Pa  at  1  meter  distance  from  the  source.  The  graphical 
results  correspond  to  the  standard  definition  of  transmission 
loss.  [Ref  7:p.  67]  The  table  in  Figure  12  shows  the  required 
input  file  blocks.  The  first  block  contains  the  title  of  the 
run,  the  second  block  the  options  selected.  For  this  study 
the  following  options  were  used  at  one  point: 

•  N  the  calculation  of  normal  stress  (equal  to  negative 
pressure  in  a  fluid) . 

•  I  Hankel  transform  integrands  are  plotted  versus 
wavenumber . 

•  T  Transmission  loss  of  the  pressure  versus  range. 

•  J  Complex  integration  contour. 

•  Z  Plot  of  the  water  column  velocity  profile. 

The  third  block  contains  the  frequency  which  was  either  7  Hz 
or  10  Hz,  and  the  Integration  contour  offset.  If  left  blank  no 
contour  offset  was  Incorporated.  A  zero  input  gave  the 
default  contour  offset  which  yields  60  dB  at  the  longest  range 
considered  in  the  FFT.  The  fourth  block  is  the  number  of 
layers  (NL)  Including  the  upper  and  lower  half spaces,  and  is 
followed  by  each  layers  upper  boundary  depth  from  l  to  NL,  the 
compressional  velocity,  the  shear  velocity,  the  compressional 
attenuation  in  dB/X,  the  shear  attenuation  in  dB/X,  (both  of 
which  are  converted  from  dB/m*kHz  to  dB/X  by  multiplying  each 
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BkKk 

PuamaUt 

Units 

Limits 

I 

TXTU  1 

UU«  of  tttB 

- 

<  SO  chat. 

II 

•ptl  op%2 

output  optioiu 

- 

<  40  char. 

III 

mq: 

•ouico  fttqueucjr 

Hs 

>  0 

com 

intcifutiou  contour  offtot 

dB/A 

COFT  >  0 

IV 

ML: 

number  of  Uyete.  incl.  helfepeecs 

- 

ML  >  2 

B: 

depth  of  intcrfece 

m 

- 

CC: 

compreanonel  speed 

m/s 

CC  >0 

CS: 

sheer  speed 

m/s 

- 

&C: 

eomptessioaal  ettenuetion 

dB/A 

1C  >0 

IS: 

sheer  ettenuetioa 

dB/A 

IS  >0 

tOi 

density 

g/cm* 

10  >  0 

16: 

rros  velue  of  interfece  tonihness 

m 

- 

CLt 

coneletion  length  of  toughness 

m 

CL>0 

V 

SO: 

source  depth  (meen  for  ertey) 

m 

- 

MS: 

number  of  sources  In  ertey 

- 

MS  >0 

OS: 

verticel  source  spedng 

m 

DS  >  0 

if: 

gtesing  eagle  of  beem 

deg 

- 

lit 

ertey  type 

- 

1  <  XI  <  5 

FO: 

focal  depth  of  beem 

m 

PO  9tSD 

VI 

101: 

depth  of  ftrst  receiver 

m 

• 

103: 

depth  of  lest  receiver 

m 

103  >1B1 

Ml: 

number  of  receivers 

- 

Ml  >  0 

Zl: 

plot  output  increment 

- 

Zl  >  0 

VII 

CMZNt 

minimum  pbese  velocity 

m/s 

CRXM  >  0 

CMZ: 

mejdmnm  phese  velocity 

m/s 

- 

VIII 

■V: 

number  of  wavenumber  samples 

MUss2" 

ICS: 

irst  oem|diag  point 

- 

ZCl  >  1 

ZC3: 

lest  sempiiag  point 

— 

ZC3  <  MV 

rlgurtt  12  SAFARI  Inputs  [Ref.  7:p.  68] 


by  its  respective  velocity  and  dividing  by  1000) ,  the  density 
in  g/cm’,  and  the  RMS  roughness  which  was  left  as  zero  for  the 
study.  For  this  study  an  average  fluid  solid  interface  depth 
of  3400  n  was  used  in  all  runs.  The  next  block  contains  the 
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source  depth  which  was  122  m  for  the  study.  The  sixth  block 
contains  the  receiver  depth  of  the  uppermost  receiver,  the 
lowermost  receiver  if  more  than  one,  the  number  of  receivers 
and  the  interim  receiver  depths.  Here  all  plots  are  of  just 
one  receiver  depth  and  have  the  required  default  place  holders 
of  1  for  the  number  of  receivers  and  interim  receiver  depths. 
The  seventh  block  contains  the  minimum  and  maximum  phase 
velocity  speed.  The  eighth  block  contains  the  number  of  FFT 
sampling  points,  the  starting  sampling  point,  and  ending 
sampling  point.  The  ninth  through  thirteenth  blocks  deal  with 
the  various  plotting  parameters.  Blocks  nine,  ten,  and 
thirteen  are  used  in  generating  the  output  data  per  the 
reference  guide  [Ref  7]  and  are  shown  in  Figure  13.  A 
representative  file  sample  using  contour  offset  is  shown  in 
Figure  14. 


Block 

PunmsUr 

Units 

Limits 

IX 

MZNi 

minimntn  tuifo  on  plots 

km 

• 

IRtli 

mudmum  tans*  on  plots 

km 

- 

kUNi 

length  of  .''•wds  fot  sU  plots 

cm 

RLEN  >  0 

EIKC: 

distancs  between  tick  marks 

km 

kZMC  >  0 

X» 

TRIM: 

minimum  transmission  loss 

dB 

TXU: 

maiimum  transmission  loss 

dB 

- 

TUM: 

length  of  vertical  TL  axes 

cm 

TLEN  >  0 

TUCi 

distance  between  tick  marks 

dB 

TZNC  >  0 

XI* 

DCVr: 

minimum  depth  for  plots 

m 

OCLO: 

maximum  depth  for  plots 

m 

- 

BCLM: 

length  of  depth  axis 

cm 

DCU  >  0 

BCZM: 

distance  between  tick  marks 

m 

DCXR  >  0 

XII* 

2RZN: 

minimum  contour  level 

dB 

mil 

maximum  contour  level 

dB 

- 

ZZRCi 

contour  level  increment 

dB 

ZZRC  >  0 

XIII* 

VLEF: 

wave  speed  at  left  border 

m/s 

nzoi 

wave  speed  at  tight  border 

m/s 

- 

niMi 

length  of  wave  speed  axis 

cm 

TIER  >  0 

rzRCt 

wave  speed  tick  mark  distance 

m/s 

TZVe  >  0 

»W: 

depth  at  upper  border 

m 

OTLO: 

depth  at  lower  border 

m 

- 

BTUi 

length  of  depth  axis 

cm 

Dru>o 

SfZNi 

dcpth'axis  tkk  mark  interval 

m 

Brn>o 

*  Oaljr  bt  opUou  A,  9,  Z  ud  T.  ’  Only  foe  options  C  snd  D.  *  Only  foi  option  C. 

*  Only  foe  option  X. 


rigur*  13  SAFARI  Plot  Parameters  [Ref.  7:p.  69] 
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SAFARI  FLOAT  0 
N  T  J 
7  0 
43 

0  0  0  0  0  0  0 

0  1544.88  -1542.97  0010 
2.0  1542.97  -1543.93  0010 

33.4  1543.93  -1542.00  0010 

36.1  1542.00  -1540.93  0010 

38.9  1540.93  -1539.59  0010 

40.2  1539.59  -1538.72  0010 

42.9  1538.72  -1537.45  0010 

52.5  1537.45  -1530.84  0010 

133.8  1530.84  -1530.14  0010 

170.2  1530.14  -1530.56  0010 

230.7  1530.56  -1528.08  0010  - 

281.6  1528.08  -1527.04  0010 

321.6  1527.04  -1526.09  0010 

406.4  1526.09  -1526.05  0010  f 

498.6  1526.05  -1523.25  0010 

599.1  1523.25  -1519.53  0010 

674.3  1519.53  -1515.17  0010 

737.4  1515.17  -1511.40  0010 

811.7  1511.40  -1503.42  0010 

908.2  1503.42  -1496.45  0010 

1046.5  1496.45  -1494.06  0010 

1165.8  1494.06  -1493.51  0010 

1403.8  1493.51  -1497.37  0010 

1806.9  1497.37  -1497.44  0010 

1811.5  1497.44  -1497.85  0010 

1835.8  1497.85  -1500.13  0010 
2041.0  1500.13  -1501.43  0010 
2176.0  1501.43  -1506.97  0010 
2611.0  1506.97  -1511.75  0010 
3018.0  1511.75  -1517.13  0010 
3400.0  1517.13  -1524.38  0010 

3405.0  1463.56  116.00  0.0771  1.5428  1.46  0 

3410.0  1472.09  122.00  0.0792  1.6568  1.47  0 

3415.0  1480.61  128.00  0.0814  1.7766  1.48  0 

3420.0  1489.11  134.00  0.0837  1.9015  1.49  0 

3430.0  1506.05  146.02  0.0881  2.1567  1.51  0 

3440.0  1522.92  158.02  0.0926  2.4256  1.53  0 

3525.0  1663.28  401.77  0.1341  0  1.69  0 

3700.0  1935.24  552.79  0.2349  0  2.00  0 

4100.0  2092.86  672.14  0.2354  0  2.18  0 

4900.0  2611.13  1074.68  0.1353  0  2.32  0  i 

6350.0  5175.0  2725.0  0.1035  0.19075  2.75  0 
122 

435.32  435.32  1  1 
50  1E8 
4096  1  4096 
0.0  25  20  10 
40  120  12  20 

Tigur*  14  SAFARI  Input  File  with  standard 
contour  offset  selected 
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V.  MODEL  ZMVE8TI6ATI0M 

In  the  SAFARI  model,  two  issues  must  be  investigated.  The 

first  is  whether  and  the  two  endpoints  of  the  sampling 

regime  alias  the  result  by  truncating  the  depth-dependent 

Greens  function  at  a  point  where  the  contribution  from  the 

integrand  is  significant.  If  the  integration  is  truncated  at 

a  point  where  the  kernel's  amplitude  is  significant  the 

discontinuity  will  give  rise  to  aliasing.  The  second  is  the 

whether  a  fine  enough  sampling  of  the  horizontal  wavenumber  k 

is  taking  place,  to  ensure  no  aliasing  is  occurring  in  the 

depth-dependent  Green's  function  (the  integrand) ,  which 

oscillates  sharply  over  the  span  of  the  horizontal 

wavenumber.  As  previously  stated,  the  contour  offset  is  used 

to  smooth  the  kernel  by  stepping  off  the  axis.  This  has  the 

effect  of  reducing  the  maximum  spatial  frequency  at  which 

aliasing  commences,  and  therefore  in  range  reduces  the 

distance  at  which  the  onset  of  aliasing  occurs. 

A.  MODEL  ALZA8ZM0 

1.  Aliasing  Background 

A  standard  definition  of  aliasing  is 

"if  a  signal  contains  no  frequency  components  above  a 
frequency  f^  the  signal  can  be  uniq[uely  represented  by 
equally  spaced  samples  if  the  sampling  frequency  f,  is 
greater  than  twice  f^.  That  is,  the  sampling  frequency 
must  satisfy  the  Inequality  f,  >  2fa„."[Ref.  12] 
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The  Inequality  is  known  as  the  Nyquist  rate.  Since  frequency 
is  analogous  to  range  and  time  analogous  to  wavenumber, 
frequency  aliasing  occurs  when  the  time  domain  is 
undersampled.  Similarly,  in  SAFARI,  range  aliasing  occurs 
when  wavenumber  space  is  undersampled.  Since  the  wavenumber 
sampling  interval  Ak»(k_„-k^)  /M,  where  M  is  the  total  number 
of  sampling  points,  either  increasing  M  or  decreasing  the 
wavenumber  span  k^u-k^  will  reduce  aliasing  in  range.  Just 
as  an  Impulse  function  in  time  domain  yields  a  constant  level 
in  the  frequency  domain,  so  to  will  sharp  peaks  in  wavenumber 
space  yield  a  constant  level  in  the  spatial  domain.  If  the 
Nyquist  rate  is  violated,  i.e.,  the  frequency  k,  is  not 
greater  than  twice  k^  the  resulting  transmission  loss  curve 
will  show  aliasing. 

2 .  Aliasing  Results 

In  SAFARI,  the  horizontal  wavenumber  k,  governs  the 
distance  at  which  the  aliasing  commences.  Given  the  rapidly 
varying  integrand  resulting  from  Eg.  3.27,  the  wavenumber 
range  must  be  large  enough  and  the  sampling  period  fine  enough 
to  allow  for  accurate  integration.  There  are  five  model  input 
parameters  governing  the  wavenumber  Integration: 

•  the  minimum  phase  velocity 

•  the  maximum  phase  velocity 

•  NW  the  number  of  wavenumber  samples,  which  must  be  an 
Integer  power  of  2,  i.e.  NW«2'*. 

•  ICl  the  first  sampling  point 
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•  IC2  the  last  sampling  point 

The  relationship  between  the  limits  of  the  horizontal 
waveniunbers  k  and  phase  speeds  C  is  k^»27rf/C^  and 
k^»2Trf/C^,  where  f  is  the  frequency  of  interest.  The 
wavenumber  interval  Ak  is  determined  by  taking  (k^-k^ )  /  M, 
where,  M  is  the  total  number  of  sampling  points.  Phrased  in 
another  fashion,  »  )c^  +  (M-l)Ak.  Then  recalling  Eg  3.26 
Ar  can  be  determined  and  the  window  range  R  «MAr  can  be 
determined.  This  is  the  maximum  range  at  which  the  integrand 
can  not  be  aliased  due  to  the  truncation.  In  this  study,  all 
graphs  start  with  ICl  equal  to  1  and  IC2  equal  to  the  last 
sampling  point  M.  To  ensure  all  compressional  and  shear 
speeds  were  incorporated,  was  chosen  well  below  the 
minimum  speed  in  the  profile.  Unless  there  are  computing 
constraints,  a  very  large  value  should  be  chosen  for  to 
ensure  truncation  does  not  occur  too  soon.  The  first  graphs 
Figure  15  and  Figure  16  show  the  transmission  loss  curve  and 
the  wavenumber  Integrand.  The  transmission  loss  curve  shows 
the  beginning  of  aliasing  just  before  30  km,  which  supports 
the  cutoff  of  the  model  at  its  predicted  value  of  29000  m. 
This  number  was  arrived  at  via  Equation  3.26.  ,i.e., 
MAr*2ir/Ak,  where  Ak  is  k^  -  k.^.  The  wavenumber  integrand 
shows  a  very  confined  area  over  which  a  very  abrupt  series  of 
oscillations  take  place  on  the  axis  of  integration.  Figure 
17  and  Figure  18  compare  the  effect  of  reducing  the  number  of 
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sanpllng  points  by  half  but  keeping  Ak  the  same  by  reducing 
k^->k^  In  half.  There  are  two  Items  to  notice.  First  the 
onset  of  aliasing  occurs  at  the  predicted  range  of  29,000  m, 
the  same  as  the  previous  result.  Second,  the  graphs  are  not 
the  same.  Indicative  of  a  reduction  In  aliasing  In  the  kernel 
when  the  higher  value  of  vas  used,  as  Is  expected  since 
the  span  In  k  has  been  reduced. 

Inspection  of  a  series  a  series  of  these  curves  allows 
a  and  which  will  provide  the  greatest  alias-free  range 
for  a  given  number  of  sampling  points.  Figure  19  where 
Ca^-740  m/s,  and  C^«1E8  m/s  (1X10*) ,  show  no  aliasing  but  the 
picture  has  converged  to  a  more  stable  version.  In  this  case, 
values  of  phase  speed  below  740  m/s  and  above  1E8  m/s  do  not 
significantly  contribute.  The  Integrand  In  Figure  20  has  been 
rescaled  automatically  by  SAFARI.  There  still  Is  a  large 
degree  of  'peakness'  and  this  will  be  Investigated  later  In 
the  contour  offset  series. 

B.  COBTOUR  OmBT 

As  previously  described  the  contour  offset  Is  designed  to 
smooth  out  the  peaks  In  the  Integrand  by  moving  away  from  the 
poles  Into  the  complex  plane.  This  Is  analogous  to  going 
along  the  side  of  the  mountain  as  opposed  to  going  over  Its 
peaks.  Since  the  Integrand  Is  composed  of  many  spatial 
frequencies  and  the  kernel  Is  smoothed  by  taking  an  offset 
contour,  this  Is  the  same  as  reducing  the  maximum  range  at 
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which  the  onset  of  aliasing  will  occur.  Given  the  number  of 
sampling  points  is  not  changed,  nor  is  the  span  of  k,  aliasing 
should  begin  at  smaller  values  of  r.  Figure  21  and  Figure  22 
show  the  onset  of  aliasing  as  moving  to  shorter  ranges  as  the 
contour  is  moved  off  the  axis.  Figures  23  through  29  show  the 
effect  of  using  the  contour  offset  for  the  43  layer  model  with 
C^>b5o  m/s  and  C^«1E8  m/s.  Figure  23  has  no  offset  and  each 
following  Figure  has  an  increase  in  the  contour  offset  by  a 
factor  of  0.25  times  the  standard  contour  offset.  With  a 
standard  offset  and  4096  points,  the  aliasing  begins  at  about 
25  km.  Doubling  the  number  of  points  to  8192  moves  the 
aliasing  out  beyond  the  limits  of  the  graph,  as  seen  in  Figure 
30  and  Figure  31. 

In  Figure  19  and  Figure  20,  was  1E8  and  was  740 
with  no  contour  offset  used.  With  a  standard  contour  offset. 
Figure  32  shows  the  change  in  the  nulls.  This  results  in  a 
small  change  in  the  transmission  loss  curve  as  shown  in  Figure 
33;  the  nulls  are  deeper  and  the  overall  shape  shows  a  few 
more  peaks.  In  the  next  series,  the  standard  contour  offset 
was  used,  and  was  held  constant  at  740  m/s,  but  was 
reduced  to  3400  m/s  (Figure  34  and  Figure  35) .  This  had  little 
effect  on  the  results  indicating  that  phase  speeds  above  3400 
m/s  contribute  little  to  the  propagation.  Figure  35  bears  this 
out  as  small  k  shows  no  significant  contribution  to 
propagation.  The  largest  effect  was  between  5  to  10  km.  The 
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next  comparison  was  in  running  the  same  profile  with  1024 
points  vice  4096.  Figure  34  shows  a  transmission  loss  curve 
almost  identical  to  the  4096  pt  curve  in  Figure  34,  with  some 
very  minor  changes  in  the  nulls  around  8  km.  The  integrand  in 
Figure  37  has  changed  dramatically.  What  were  originally 
spikes  in  Figure  20,  now  have  some  width  and  continuity, 
making  for  an  easier  evaluation  of  the  integrand  and  thus  less 
aliasing  due  to  the  kernel.  Next  is  the  comparison  of  Figure 
36  with  the  8192  point  graph  in  Figure  31  both  using  the 
contour  offset  feature,  to  observe  the  losses  due  to  the 
change  in  parameters.  Clearly  not  much  is  lost  in  the  overall 
shape  of  the  transmission  loss  curve.  A  few  more  features  are 
present  with  deeper  nulls  and  smoother  peaks  and  some 
resolution  of  a  few  peaks  around  6  km,  but  otherwise  the  shape 
is  the  same,  with  the  1024  point  graph  taking  less  time  to 
produce.  Finally  increasing  the  number  of  sampling  points  to 
16384  points  (Figure  38)  yielded  no  change  in  the  graph  over 
that  of  the  8192  point  graph  (Figure  31)  over  the  range  of 
Interest.  At  this  point  the  model  has  converged.  Note  the 
difference  in  the  Integrands  between  Figure  39  and  Figure  37. 

Finally,  Figure  40  shows  the  difference  between  the  water 
column  alone  and  using  the  water  column  and  sediment  combined. 
Clearly  sediment  contribution  to  the  final  transmission  loss 
is  important. 
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VX.  COMPARISON  OF  SAFARI  WITH  SWALLOW  FLOAT  0  DATA 

In  this  chapter  the  results  of  the  Swallow  Float  records 
are  presented  graphically  for  7  Hz  and  10  Hz  both  without  and 
with  range  offsets  from  regions  1  and  2.  The  results  are 
visually  compared  and  contrasted.  Next  the  data  is  compared 
with  the  model  results.  In  all  cases  Matlab  386  was  used  for 
the  graphic  presentations.  All  figures  are  located  in 
Appendix  0. 

A.  7  HZ  DATA 

Figure  41  shows  the  7  Hz  transmission  loss  data  taken 
during  region  1.  Figure  42  shows  the  7  Hz  transmission  loss 
data  taken  during  region  2.  Figure  43  compares  region  l  and 
region  2  data  by  offsetting  region  1  to  match  predominant 
features.  This  offset  was  2.4  km. 

B.  10  HZ  DATA 

Figure  44  shows  the  10  Hz  transmission  loss  data  taken 
from  region  1.  Figure  45  shows  the  10  Hz  transmission  loss 
data  taken  from  region  2.  Figure  46  compares  region  l  with 
region  2  by  offsetting  region  l  to  match  predominant  features. 
This  offset  was  2.5  km. 

C.  SAFARI  MODEL  VS  7  HZ  DATA 

Figure  47  shows  the  SAFARI  model  at  7  Hz  vs  7  Hz 
transmission  loss  data  taken  from  region  1.  Figure  48  shows 
the  SAFARI  model  at  7  Hz  vs  7  Hz  transmission  loss  data  taken 
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from  region  2.  Figures  49  and  50  compare  the  SAFARI  model  at 
7  Hz  vs  7  Hz  transmission  loss  data  taken  from  regions  1  and 
2  respectively,  offset  in  range  to  match  predominant  features. 
The  offset  was  2.4  km  for  region  1  and  none  for  region  2. 

D.  SAFARI  MODEL  VS  10  B2  DATA 

Figure  51  shows  the  SAFARI  model  at  10  Hz  vs  lO  Hz 
transmission  loss  data  taken  from  region  1.  Figure  52  shows 
the  SAFARI  model  at  10  Hz  vs  10  Hz  transmission  loss  data 
taken  from  region  2.  Figures  53  and  54  compare  the  SAFARI 
model  at  10  Hz  vs  10  Hz  transmission  loss  data  taken  from 
regions  1  and  2  respectively,  offset  in  range  to  match 
predominant  features.  The  offset  was  0.5  km  for  region  1  and  - 
0.75  )cm  for  region  2. 
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VII.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

The  Initial  findings  in  the  Swallow  float  0  data  show 
excellent  agreement  with  the  SAFARI  model;  the  model  produces 
*  a  result  that  compares  well  with  ocean  data.  The  model  shows 

,  that  there  is  a  distinct  problem  with  the  range  used  for  the 

data  recorded  in  region  l.  No  explanation  was  found  as  to  why 
the  offset  found  in  region  1  at  7  Hz  was  much  larger  than  the 
navigation  error  mentioned  in  Reference  4.  The  data  at  7  Hz 
in  region  2  and  at  10  Hz  in  both  regions  l  and  2  match  up 
reasonably  well  given  errors  in  the  navigation.  The  slope  of 
the  bottom  appears  to  have  very  little  effect  on  the 
transmission  loss  in  this  case.  By  comparing  results  obtained 
from  modelling  the  velocity  profile  of  the  water  column  alone 
versus  one  with  sediments  included,  the  impact  of  the  bottom 
and  the  importance  of  modelling  the  bottom  become  apparent. 

SAFARI  is  extremely  sensitive  to  aliasing  based  on  its 
i>  input  values  of  phase  speed,  number  of  sampling  points,  and 

limits  of  integration.  The  contour  offset  does  provide  a  way 
to  generate  less  aliasing  from  poles  at  the  cost  of  aliasing 
at  long  ranges.  This  can  be  corrected  by  increasing  the 
number  of  sampling  points. 
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B.  RECOMMENDATIONS 

Two  areas  which  can  be  addressed  in  future  theses  are 
investigation  of  data  obtained  on  other  floats  and 
investigation  of  the  effect  the  starting  bottom  depth  has  on 
the  model. 
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APPENDIX  A 


A.  NAVIGATIONAL  DATA 

The  record  of  the  position  of  the  R/V  Gyre  is  enclosed. 
It  is  uncorrected  with  respect  to  the  navigational  offset 
mentioned  in  Chapter  2.  The  last  page  has  the  time  the 
changes  in  navigational  source  occurred. 
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30.9353 

-74.95 

4.525 

15.22939 

8/05/90  07K)SKX) 

30.9411 

-74.9483 

4.3338 

15.91076 

8/05/90  07:10K)0 

30.9467 

-74.9469 

4.09857 

14.03162 

8/05/90  07:15.-00 

30.9522 

-74.9456 

4.09911 

14.03162 

8/05/90  07:20.-00 

30.9578 

-74.9447 

4.03628 

8.51249 

8/05/90  07:25.-00 

30.9631 

-74.9439 

3.83785 

8.95455 

8/05/90  07:30:00 

30.9683 

-74.9431 

3.83874 

9.03535 

8/05/90  07:35K)0 

30.9742 

-74.9408 

4.44119 

20.83887 

8/05/90  07:40.-00 

30.98 

-74.9389 

4.3886 

18.44057 

8/05/90  07:45.00 

30.9858 

-74.9369 

4.38828 

18.44619 

30.9914 

-74.9339 

4.47225 

28.78662 

8/05/90  07:55.-00 

30.9972 

-74.9314 

4.5119 

23.21441 

8/05/90  080000 

31.0033 

-74.9292 

4.6389 

20.02904 

8/05/90  080500 

31.0092 

-74.9272 

4.39068 

18.37872 

8/05/90  08:10.-00 

31.0153 

-74.925 

4.6413 

20.02904 

8/05/90  08:1500 

31.0211 

-74.9228 

4.4512 

20.83887 

31.0369 

-74.9283 

4.59203 

25.46043 

■  V'tVXv  'tvM 

31.0428 

-74.9256 

4.59503 

25.45316 

■  I'V.nW 

31.0486 

-74.9228 

4.59548 

25.46043 

1  8/05/90  08:4500 

31.0542 

-74.92 

4.41593 

26.55718 
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8/0S/90  09:10K)0 
8/0S/90  09:1SK)0 
8/05/90  09:20.*00 
8/05/90  09:25.-00 
8/05/90  09:30K)0 
8/05/90  09:35K)0 
8/05/90  09:40K)0 
8/05/90  09:45K)0 
8/05/90  09:50.‘00 
8/05/90  09:55.-00 
8/0S/9O  lOrOOKX) 
8/05/90  10K)5.‘00 
8/05/90  10:10.-00 
8/05/90  10:15K)0 
8/05/90  10:20.-00 
8/05/90  10:25.-00 
8/05/90  10:30:00 
8/05/90  10:35KX) 
8/05/90  10:50.-00 
8/05/90  10:55.-00 
8/05/90  IIKXIKK) 
8/05/90  11K)5K)0 
8/05/90  11:10K)0 
8/05/90  ll:15.-00 
8/05/90  ll-JOKX) 
8/0S/9O  11:25.-00 
8/05/90  11:30:00 
8/05/90  11:35K)0 
8/05/90  11:40KX) 
8/05/90  11:45K)0 
8/05/90  11:50K)0 
8/05/90  11:55:00 
8/05/90  12:00K)0 
8/05/90  12.-05K)0 
8/05/90  12:10K)0 
8/05/90  12:15K)0 
8/05/90  12:20.-00 
8/05/90  12:25K)0 
8/05/90  12:30K)0 
8/05/90  12:35KX) 
8/05/90  12:40.-00 
8/05/90  12:45.-00 
8/05/90  12:50K» 
8/05/90  12:55:00 
8/05/90  1300K)0 
8/05/90  13K)5K)0 
8/05/90  13:10:00 


31.0842 

31.0903 

31.0961 

31.1022 

31.1086 

31.1147 

31.1211 

31.1275 

31.1317 

31.1364 

31.1422 

31.145 

31.1481 

31.1508 

31.1547 

31.1583 

31.1622 

31.1661 

31.1942 

31.1994 

31.2047 

31.21 

31.2153 

31.2206 

31.2258 

31.2314 

31.2369 

31.2425 

31.2481 

31.2539 

31.2597 

31.2656 

31.2711 

31.2769 

31.2825 

31.2881 

31.2939 

31.2994 

31.3053 

31.3111 

31.3169 

31.3228 

31.3289 

31.3347 

31.3406 

31.3464 

31.3522 


-74.9089 

-74.9069 

-74.905 

-74.9042 

-74.9031 

-74.9022 

-74.9014 

-74.9006 

-74.9106 

-74.9194 

-74.9242 

-74.9197 

-74.9153 

-74.9106 

-74.9092 

-74.9078 

-74.9064 

-74.9047 

-74.8933 

-74.8914 

-74.8892 

-74.8872 

-74.8842 

•74.8814 

-74.8786 

-74.8764 

•74.875 

-74.8739 

-74.8728 

-74.8714 

-74.8692 

-74.8667 

-74.8647 

-74.8631 

•74.8614 

-74.8597 

-74.8581 

-74.8567 

-74.8553 

-74.8547 

-74.8539 

-74.8531 

•74.8522 

-74.8514 

-74.8506 

-74.8497 

-74.8492 


4.4368 

4.43629 

4.59581 

4.40502 

4.43673 

4.6636 

4.43702 

4.6347 

4.6362 

7.53516 

7.02376 

5.32476 

3.67764 

3.78748 

3.83816 

2.96377 

2.77359 

2.96266 

3.03103 

4.10741 

4.03866 

4.11006 

4.03971 

4.3702 

4.27381 

4.27721 

4.29993 

4.11826 

4.0766 

4.0778 

4.31387 

4.48658 

4.563 

4.23243 

4.36657 

4.17143 

4.17455 

4.36449 

4.12179 

4.31547 

4.21842 

4.2443 

4.2424 

4.44027 

4.24171 

4.24311 

4.24179 

4.21855 


7.8 1942 

7.74921 

17.65898 

18.44619 

7.74921 

9.8889 

7.74921 

7.41751 

7.41532 

292.6106 

298.0035 

321.0036 

58.00346 

55.46642 

59.50205 

19.74811 

21.03548 

19.64848 

23.15449 

22.81514 

20.23541 

22.80774 

20.23541 

30.1004 

27.68852 

27.75342 

21.84887 

14.03624 

11.26075 

11.33641 

13.39081 

20.83887 

23.21441 

19.30407 

15.98002 

16.67036 

16.73713 

15.9157 

14.03162 

13.39081 

5.45448 

8.10913 

8.1878 

7.74921 

8.11436 

8.11174 

8.11436 

5.45448 
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■  IP  \mm 

1  miinM 

8/05/90  13:35K)0 
8/05/90  13:40:00 
8/05/90  13:50:00 
8/05/90  13:55:00 
8/05/90  14:00:00 
8/05/90  14:15:00 

31.3814 

31.3875 

31.3986 

31.4047 

31.4106 

31.4261 

•74.8486 

-74.8478 

-74.8644 

-74.8636 

-74.8628 

-74.8586 

3.82097 

4.44054 

4.20444 

4.44057 

4.24199 

2.60806 

6.0241 

7.74921 

2.non 

7.74921 

8.11436 

32.48508 

8/04/90 . 

,.1615-1720 

LORAN 

8/04-05/90  . 

.1725-0725 

GPS 

8/05/90 . 

..0730-1025 

SAT  NAV 

8/05/90 . 

..1030-1320 

LORAN 

8/05/90 . 

..1325-1415 

GPS 
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APPEHDIX  B 


A.  SOUND  VELOCITY  PROFILES 

Representative  samples  of  the  sound  velocity  profile  of 
the  water  column  used  in  the  SAFARI  model  are  presented  in 
this  appendix.  A  full  listing  of  all  the  samples  made  is 
available  in  Reference  1. 
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SITE  **0**  PROVINCE  PROFILES: 

SITE  -D-  PROVINCE  ONE  SITE  "D"  PROVINCE  TWO 

31.4833  «75J644  90  8  3850  40  1  3US91  -74.9573  90  g  3gS0  32  2 


TERUimiP 
(M)  Ideoc 

TSEPTOTSaP" 

(M)  DEGC 

■1  \^\mm 

■I  \v\wm 

26.3 

27.52 

36.84 

1542.99 

10 

28.00 

36.63 

1543.40 

29.1 

26.96 

36.84 

1541.80 

15.7 

28.00 

36.77 

1543.78 

43.6 

26.16 

37.44 

1540.91 

311 

27.30 

37.44 

1543.26 

51.7 

25.32 

37.43 

1539.12 

34.8 

26.90 

37.44 

1542.42 

54.1 

24.94 

37.42 

1538.28 

37.5 

26.00 

37.44 

1540.46 

54.7 

24.69 

37.39 

1537.68 

57.9 

23.90 

37.24 

1535.67 

59.8 

24.35 

37.16 

1536.68 

81.0 

22.90 

36.78 

1533.09 

67.6 

24.12 

36.90 

1535.97 

87.8 

22.30 

36.75 

1531.66 

80.9 

23.24 

36.78 

1533.94 

108.1 

21.80 

36.75 

1530.72 

9S.9 

22.45 

36.75 

153117 

128.4 

20.90 

36.72 

1528.67 

99.7 

22.42 

36.76 

153116 

205.2 

19.50 

36.70 

1526.13 

104J 

22.01 

36.76 

1531.21 

418.3 

18.00 

36.47 

1525.13 

114.6 

21.73 

36.74 

1530.64 

480.2 

17.30 

36.41 

1524.04 

118.3 

21.40 

36.74 

1529.83 

562.7 

15.70 

36.22 

1520.33 

134.3 

21.10 

36.71 

1529.28 

630.3 

14.80 

35.98 

1518.33 

162.3 

19.97 

36.71 

1526.71 

685.9 

13.60 

35.80 

1515.14 

206.3 

19.26 

36.70 

1525.50 

709.1 

12.80 

35.73 

1512.78 

267J 

19.00 

36.61 

1525.66 

789.9 

10.80 

35.53 

1506.95 

428.0 

17.07 

36.46 

152155 

812.9 

10.40 

35.45 

1505.79 

561.8 

14J5 

36.22 

1516.06 

965.0 

7.60 

35.12 

1497.47 

622.0 

13.57 

36.02 

1514.26 

1091.4 

6.00 

35.08 

1493.23 

668.1 

12.18 

35.86 

1510.16 

1238.4 

5.20 

35.05 

1492.42 

691.3 

11.86 

35.78 

1509.35 

1478.4 

4.50 

35.01 

1493.51 

763.8 

9.99 

35.60 

1503.68 

1789.5 

4.20 

34.98 

1497.49 

868X) 

8.02 

35.29 

1497.69 

1835.8 

4.00 

34.98 

1497.44 

900.0 

7.75 

35  J3 

1497.10 

2041.0 

3.81 

34.98 

1500.13 

913.0 

7.41 

35.19 

1495.98 

2176.0 

3.58 

34.96 

1501.43 

989.0 

6.50 

35.10 

1493.54 

2611.0 

3.13 

34.95 

1506.97 

1090.0 

6.04 

35.08 

1493.37 

3018.0 

2.62 

34.92 

1511.75 

1160.0 

5.64 

35.07 

149193 

3400.0 

2J3 

34.90 

1517.13 

129a0 

5.02 

35.04 

149154 

3850.0 

2.18 

34.88 

1524.38 

1416.0 

4.60 

35.01 

149188 

1628.0 

4.19 

34.99 

1494.70 

2041.0 

3.81 

34.98 

isoai3 

2176.0 

3J8 

34.96 

1501.43 

2611.0 

3.13 

34.95 

1506.97 

3018.0 

2.62 

34.92 

1511.75 

34oao 

133 

34.90 

1517.13 

3850.0 

118 

34.88 

1524.38 
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SITE  *0"  PROVINCE  THREE  SITE  "D"  PROVINCE  FOUR 
3a.7<15  -75.0077  90  t  3850  32  3  31.1660  -74.9048  90  8  3850  32  4 


KSii 

TSEPramiri 

(M)  DEGCl 

lllfITiM 

■1 

2.0 

27.80 

36.63 

1342.97 

2.0 

27.80 

36.63 

1542.97 

19.8 

27.90 

36.81 

1343.68 

■33.4 

27.60 

37.44 

1343.93 

32.1 

27.40 

37.44 

1343.48 

36.1 

26.70 

37.44 

1342.00 

36.1 

27.00 

37.44 

1342.67 

38.9 

26.20 

37.44 

1340.93 

38.9 

26.40 

37.44 

1541.38 

40.2 

25.60 

37.44 

1539.59 

47.0 

23.30 

37.44 

1339.02 

42.9 

23.20 

37.44 

1338.72 

67.4 

24.30 

36.90 

1336.41 

32.5 

24.60 

37.42 

1537.45 

\77J0 

22.20 

36.72 

1332.02 

'133.8 

21.70 

36.71 

1330.84 

2173 

20.00 

36.69 

1327.68 

170.2 

21.20 

36.71 

1330.14 

2S4.8 

19.40 

36.63 

1326.38 

230.7 

21.00 

36.67 

1330.56 

AXnA 

18.40 

36.47 

1326.11 

281.6 

19.80 

36.39 

1328.08 

347.0 

17.00 

36.27 

1324.08 

321.6 

19.20 

36.33 

1327.04 

649.7 

13.20 

33.92 

1319.84 

406.4 

18.40 

36.48 

1526.09 

667.8 

14.30 

33.86 

1317.84 

498.6 

17.90 

36.37 

1326.05 

106^ 

13.30 

33.73 

1313.08 

599.1 

16.50 

36.11 

1323.25 

786.1 

11.90 

33.34 

1310.73 

674.3 

13.00 

33.83 

1319.33 

839.7 

11.20 

35.35 

1508.95 

737.4 

13.40 

35.67 

1515.17 

9313 

8.80 

33.14 

1501.82 

811.7 

12.00 

35.43 

1311.40 

988.9 

7.80 

35.10 

1498.61 

908.2 

9.40 

33.21 

1303.42 

1037.8 

7,10 

33.09 

1496.70 

1046J 

7.00 

33.09 

1496.43 

1136.1 

6.00 

33.07 

1493.97 

1163.8 

3.90 

35.07 

1494.06 

12993 

3.10 

33.04 

1493.02 

-1403.8 

4.80 

33.02 

1493.31 

1498.7 

4.30 

35.00 

1493.85 

1806.9 

4.10 

34.98 

1497.37 

1783.7 

4.10 

34.98 

1496.97 

1811J 

4.10 

34.98 

1497.44 

1833.8 

4.10 

34.98 

1497.83 

1833.8 

4.10 

34.98 

1497.83 

2041.0 

3.81 

34.98 

1300.13 

2041i> 

3.81 

34.98 

1300.13 

2176.0 

3.38 

34.96 

1301.43 

2176.0 

3.38 

34.96 

1501.43 

2611.0 

3.13 

34.93 

1306.97 

2611.0 

3.13 

34.95 

1306.97 

3018.0 

2.62 

34.92 

1311.73 

3018.0 

2.62 

34.92 

1311.73 

3400.0 

2.33 

34.90 

1517,13 

3400.0 

2.33 

34.90 

1517.13 

3830.0 

2.18 

34.88 

1324.38 

3830.0 

2.18 

34.88 

1324.38 
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(133J)  Mld30  HU» 


DEPTH  (MnCRS) 


mpth  (Mf 


depth  (ren) 


APPENDIX  C 


A.  6E0AC0U8TIC  BOTTOM 

The  geoacoustic  bottom  model  and  its  results  are  presented 
in  this  appendix  and  are  taken  from  Reference  l.  The 
conversion  of  the  compressions!  and  shear  attenuation 
coefficients  from  dB/(m*kHz)  to  dB/A,  needed  as  an  input  in 
SAFARI,  is  as  follows: 

A.=A„*C,/1000 

where  A,  is  the  attenuation  coefficient  of  the  respective 
compressions!  wave  or  shear  wave  in  the  SAFARI  model,  is 
the  attenuation  coefficient  of  the  geoacoustic  model  C,  is  the 
velocity  of  the  x  wave,  where  x  is  compressions!  or  shear,  and 
1000  is  a  scaling  factor  from  kHz  to  Hz. 
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Unconsolidated  layer  •  Calcareous  Sediment,  (DSDP  Well  533) 

Vp  •  1.455  +  1.713(D)  -  0.374(0)2  (Hamilton.  1980) 

Vs  «  110.0  *  1.2005(D)  (for  Vp  <  1.512)  (Hamilton,  1976) 

Vs  « (3.884  X  Vp)  •  5.757  (for  Vp  1.512-1.555  km/s)  (Hamilton.  1980) 
Vs  *  0.991  -  (1.136  X  Vp)  +  0.47(Vp)2  (for  Vp  1.65-2.15  km/s) 

Vs  •  (0.78  xVp)- 0.962  (for  Vp>  2.15  km/s) 

Kp:  Estimated  from  Hamilton,  1980 
Ks:  Estimated  from  Hamilton,  1980 

p  -  1.1379  (Vp)  -  0.20484  (Hamilton.  1978) 

Gas  layer  beneath  BSR  (m/s)- 

Vp  «  1090.4  +  1.2892(D)  -  [(2.8628  xlO^  )x(D^] 

Vs  »  Vp/1.9  (Hamilton,  1980;  Pandit  and  King,  1982) 

Kp:  Estimated  from  Hamilton,  1980 
Ks:  Estimated  from  Hamilton,  1980 

p  «  0.1 12  (Vp)  *  2.088  (HamUton,  1978) 


BWSS  «  Bottom  Water  Sound  Speed 
D<*Depdi 

Vp  «  Coo^fessional  Velocity 
Vs  ■  Shear  Velocity 
p«Density 

I4>  *  Comptessional  Attenuation 
Ks  ■  Shear  Attenuadoo 


63 


Geoacoustic  models  for  Site  ”D" 

Event  lUg  I  Start  31*26*13"  N.  74*49'27'*  W 
Bearing  193*  Range  0  nm 


D 

(m) 

(mfs) 

- Ts - 

(m/s) 

- Kp - 

(dB/m*kHz) 

— rs — 

(dB/m*lcHr) 

mm 

U.U 

3.0 

1463.36 

116.00 

0.0527 

13.30 

1.46 

10.0 

1472.09 

122.00 

0.0538 

13.58 

1.47 

15.0 

1480.61 

128.00 

0.0550 

13.88 

1.48 

20.0 

1489.11 

134.00 

0.0562 

14.19 

1.49 

30.0 

1306.03 

146.02. 

0.0585 

14.77 

1.31 

40.0 

1522.92 

158.02 

0.0608 

15.35 

1.53 

30.0 

1339.72 

223.27 

0.0632 

15.93 

1.35 

73.0 

1581.37 

313.02 

0.0690 

17.42 

1.39 

100.0 

1622.56 

359.85 

0.0748 

18.88 

1.64 

125.0 

1663.28 

401.77 

0.0806 

20.00 

1.69 

130.0 

1703.34 

419.74 

0.0865 

20.00 

1.73 

200.0 

1782.64 

459.49 

0.0981 

20.00 

1.82 

230.0 

1859.88 

503.98 

0.1098 

20.00 

1.91 

300.0 

1933.24 

552.79 

0.1214 

20.00 

2.00 

400.0 

2080.36 

661.82 

0.1575 

20.00 

2.10 

500.0 

2218.00 

768.04 

0.1425 

20.00 

2.17 

600.0 

2348.16 

869.36 

0.1275 

20.00 

2.21 

630.0 

1807.41 

458.71 

0.2400 

30.00 

1.67 

673.0 

1951.94 

366.07 

0.1762 

25.00 

1.92 

700.0 

2092.86 

672.14 

0.1125 

20.00 

2.18 

800.0 

2133.34 

719.32 

0.0975 

20.00 

2.25 

900.0 

2213.74 

766.28 

0.0825 

20.00 

2.26 

1000.0 

2273.00 

8 12  JO 

0.0675 

17.04 

2.27 

1300.0 

2611.13 

1074.68 

0.0518 

13.08 

2.32 

2000.0 

3153.00 

1498.90 

0.0361 

9.11 

2.37 

2430.0 

3836.13 

2030.18 

0.0220 

5.55 

2.42 

2300.0 

3861.13 

2049.68 

0.0205 

5.17 

2.43 

2675.0 

3948.63 

2117.93 

0.0150 

3.79 

2.44 

2700.0 

3961.13 

2127.68 

0.0142 

3.58 

2.43 

2800.0 

4011.13 

2166.68 

0.0111 

2.80 

2.46 

Basement; 

Vp:  S17S  m/t 

Vs:  2723  m^s 

Ks:  0.07  dBAn«kHz 

fc  2.73  gfcc 


Appendix  G.  Event  lUg  IStart  31^6  '13**  N.  1V*A9  *27**  W 
Bearing  193*  Range  10.45  nm 


(<lB/m*kHz)  (dB/m*kHz) 


ifcftm 


Basement: 
Vk  51* 
27; 

Kp:  0.C 
Ks:  0.C 
p:  2.7 


5175  m/s 
2725  m/s 
0.02dB/in*kHz 
0.07dB/lBi*kHz 
2.75  g/oe 


13.30 

1.46 

13.58 

1.47 

13.88 

1.48 

14.19 

1.49 

14.77 

1.51 

15.35 

1.53 

15.95 

1.55 

17.42 

1.59 

18.88 

1.64 

20.00 

1.69 

20.00 

1.73 

20.00 

1.82 

20.00 

1.91 

20.00 

2.00 

20.00 

2.10 

20.00 

2.17 

30.00 

1.66 

25.00 

1.90 

20.00 

2.14 

20.00 

2.18 

20.00 

2.25 

20.00 

2.26 

17.04 

2.27 

1  3.08 

2.32 

9.11 

2.37 

5.55 

2.42 
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Appendix  G.  Event  ILeg  IStart  31*26  *13"  N.  74*49  ’ll”  W 
Bearing  193*  Range  1421  nm 


- D - 

(m) 

- Tv - 

(m/$) 

- Vs - 

(tn/s) 

- Kp - 

(dB/m*kHz) 

- Ks - 

(dB/in*kHz) 

P 

g^cc 

0.0 

r.4j 

5.0 

1463.56 

116.00 

0.0527 

13.30 

1.46 

10.0 

1472.09 

122.00 

0.0538 

13.58 

1.47 

15.0 

1480.61 

128.00 

0.0550 

13.88 

1.48 

20.0 

1489.11 

134.00 

0.0562 

14.19 

1.49 

30.0 

1506.05 

146.02. 

0.0585 

14.77 

1.51 

40.0 

1522.92 

158.02 

0.0608 

15.35 

1.53 

50.0 

1539.72 

223.27 

0.0632 

15.95 

1.55 

75.0 

1581.37 

313.02 

0.0690 

17.42 

1.59 

100.0 

1622.56 

359.85 

0.0748 

18.88 

1.64 

125.0 

1663.28 

401.77 

0.0806 

20.00 

1.69 

150.0 

1703.54 

419.74 

0.0865 

20.00 

1.73 

200.0 

1782.64 

459.49 

0.0981 

20.00 

1.82 

250.0 

1859.88 

503.98 

0.1098 

20.00 

1.91 

300.0 

1935.24 

552.79 

0.1214 

20.00 

2.00 

400.0 

1560.27 

330.91 

0.3150 

30.00 

1.58 

425.0 

1659.06 

389.67 

0.2120 

25.00 

1.69 

500.0 

1953.88 

565.68 

0.1425 

20.00 

2.02 

580.0 

2012.19 

608.14 

0.1305 

20.00 

2.08 

600.0 

2026.15 

618.78 

0.1275 

20.00 

2.10 

650.0 

2060.11 

645.42 

0.1200 

20.00 

2.14 

700.0 

2092.86 

672.14 

0.1125 

20.00 

2.18 

800.0 

2155.54 

719.32 

0.0975 

20.00 

2.25 

900.0 

2215.74 

766.28 

0.0825 

20.00 

2.26 

1000.0 

2275.00 

812.50 

0.0675 

17.04 

2.27 

1500.0 

2611.13 

1074.68 

0.0518 

13.08 

2.32 

2000.0 

3155.00 

1498.90 

0.0361 

9.11 

2.37 

2450.0 

3836.13 

2030.18 

0.0220 

5.55 

2.42 

Basement: 

Vp;  5175  m/i 

Vs:  2725  m/s 

Kp:  0.02dB/m*kHz 
Ks:  0.07  dBAn*kHz 

p:  2.75  g/oc 
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Appendix  G.  Event  lUg  IStait  31*26  *13"  N.  74*49  77"  W 
Bearing  193*  Range  28.63  nm 


- E5 - 

(m) 

(mS) 

- Vs - 

(m/s) 

- K5 - 

(dB/m*kHz) 

9 

t/ce 

. . Ilf 

Hi 

5.0 

1463.56 

116.00 

0.0527 

BeH 

10.0 

1472.09 

122.00 

0.0538 

15.0 

1480.61 

128.00 

O.OSSO 

20.0 

1489.11 

134.00 

0.0562 

14.19 

1.49 

30.0 

1506.05 

146.02. 

0.0585 

14.77 

1.51 

40.0 

1522.92 

158.02 

0.0608 

15.35 

1.53 

50.0 

1539.72 

223.27 

0.0632 

15.95 

1.55 

75.0 

1581.37 

313.02 

0.0690 

17.42 

1.59 

100.0 

1622.56 

359.85 

0.0748 

18.88 

1.64 

125.0 

1663.28 

401.77 

0.0806 

20.00 

1.69 

150.0 

1703.54 

419.74 

0.0865 

20.00 

1.73 

200.0 

1782.64 

459.49 

0.0981 

20.00 

1.82 

250.0 

1859.88 

503.98 

0.1098 

20.00 

1.91 

300.0 

1935.24 

552.79 

0.1214 

400.0 

2080.36 

661.82 

0.1575 

20.00 

500.0 

2218.00 

768.04 

0.1425 

20.00 

2.17 

600.0 

2348.16 

869 J6 

0.1275 

20.00 

2.21 

700.0 

2470.84 

965.26 

0.1125 

20.00 

2.24 

800.0 

2586.04 

1055.11 

0.0975 

20.00 

2.26 

900.0 

2693.76 

1139.13 

0.0825 

20.00 

2.28 

1000.0 

2794.00 

1217.32 

0.0675 

17.04 

2.29 

13.08 

2.34 

KtruTiS 

1678J0 

9.11 

2.36 

3955.00 

2122.90 

5.15 

2.42 

4385.00 

2458.30 

1.77 

2.47 

Basement: 

Vp:  5175  B^s 

Vs:  2725  nVs 

Kp:  0.02dBAn*kHz 
Ks:  0.07dB/in*Ulz 
p:  2.75  g/bc 
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Appendix  O.  Event  ILeg  IStait  31^6  *13"  N,  74*'49  77"  W 
Beaiing  193**Range  40.13  nm 


D - 

(m) 

Vp 

(m/s) 

— Vi — 

(m/%) 

- Kp - 

(dB/ni*kHz) 

Ks 

(<lB/m*kHzl 

m 

■  1  Ilf 

W  iTilHI 

5.0 

1463.56 

116.00 

0.0527 

13.30 

1.46 

10.0 

1472.09 

122.00 

0.0538 

13.58 

1.47 

15.0 

1480.61 

128.00 

0.0550 

13.88 

1.48 

20.0 

1489.11 

134.00 

0.0562 

14.19 

1.49 

30.0 

1506.05 

146.02. 

0.0585 

14.77 

1.51 

40.0 

1522.92 

158.02 

0.0608 

15.35 

1.53 

30.0 

1339.72 

223.27 

0.0632 

15.95 

1.55 

75.0 

1581.37 

313.02 

a0690 

17.42 

1.59 

100.0 

1622.56 

359.85 

0.0748 

18.88 

1.64 

125.0 

1663.28 

401.77 

0.0806 

20.00 

1.69 

150.0 

1703.54 

419.74 

0.0865 

20.00 

1.73 

200.0 

1782.64 

459.49 

0.0981 

20.00 

1.82 

250.0 

1859.88 

503.98 

0.1098 

20.00 

1.91 

300.0 

1935.24 

552.79 

0.1214 

20.00 

2.00 

400.0 

2080.36 

661.82 

0.1575 

20.00 

2.10 

500.0 

2218.00 

768.04 

0.1425 

20.00 

2.17 

580.0 

1741.37 

424.64 

0.2610 

30.00 

1.65 

605.0 

1952.21 

503.44 

0.2110 

25.00 

1.82 

650.0 

2060.11 

645.42 

0.1200 

20.00 

2.14 

700.0 

2092.86 

672.14 

0.1125 

20.00 

2.18 

800.0 

2155.54 

719.32 

0.0975 

20.00 

2.25 

900.0 

2215.74 

766.28 

0.0825 

20.00 

2.26 

1000.0 

2273.00 

812.50 

0.0675 

17.04 

2.27 

1500.0 

2611.13 

1074.68 

0.0518 

13.08 

2.32 

2000.0 

3135.00 

1498.90 

0.0361 

9.11 

2.37 

2450.0 

3836.13 

2030.18 

0.0220 

5.55 

2.42 

2500.0 

3861.13 

2049.68 

0.0205 

5.17 

2.43 

2675.0 

3948.63 

2117.93 

0.0150 

3.79 

2.44 

2700.0 

3961.13 

2127.68 

0.0142 

3.58 

2.45 

2800.0 

4011.13 

2166.68 

0.0111 

2.80 

2.46 

3000.0 

4111.13 

2244.68 

0.0048 

1.21 

2.48 

Basemenc 
Vp:  317S  ta/i 

Vs:  272S  xal% 

Kp:  0.02dBAn*kHz 
Ks:  0.07dB/iai«kHz 
p:  2.73  f/cc 
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Appendix  G.  Event  ILeg  ISun  31^6  *13"  N.  74*49  *27"  W 
Bearing  193*  Range  42  nin 


mm 

mm 

- Kp - 

((lBAtn*kHz) 

- Ki - 

(dB/in*kHz) 

1463.36 

10.0 

1472.09 

13.0 

1480.61 

128.00 

0.0530 

13.88 

1.48 

20.0 

1489.11 

134.00 

0.0362 

14.19 

1.49 

30.0 

1306.03 

146.02. 

0.0585 

14.77 

1.51 

40.0 

1322.92 

138.02 

0.0608 

15.35 

1.53 

30.0 

1339.72 

223.27 

0.0632 

13.95 

1.53 

73.0 

1381.37 

313.02 

0.0690 

17.42 

1.59 

100.0 

1622.36 

339.83 

0.0748 

18.88 

1.64 

123.0 

1663.28 

401.77 

0.0806 

20.00 

1.69 

130.0 

1703.34 

419.74 

0.0865 

20.00 

1.73 

200.0 

1782.64 

439.49 

0.0981 

20.00 

1.82 

230.0 

1839.88 

303.98 

0.1098 

20.00 

1.91 

300.0 

1933.24 

352.79 

0.1214 

20.00 

2.00 

400.0 

2080.36 

661.82 

0.1575 

20.00 

2.10 

300.0 

1663.30 

384.02 

0.2830 

30.00 

1.63 

323.0 

1734.48 

442.78 

0.2436 

23.00 

1.75 

600.0 

2026.13 

618.78 

0.1273 

20.00 

2.10 

700.0 

2092.86 

672.14 

0.1123 

20.00 

2.18 

800.0 

2133.34 

719.32 

0.0973 

20.00 

2.25 

900.0 

2213.74 

766.28 

0.0823 

20.00 

2.26 

1000.0 

2273.00 

812.30 

0.0673 

17.04 

2.27 

1300.0 

2611.13 

1074.68 

0.0318 

13.08 

2.32 

2000.0 

3135.00 

1498.90 

0.0361 

9.11 

2.37 

2300.0 

3861.13 

2049.68 

0.0203 

3.i7 

2.43 

2117.93 

0.0130 

3.79 

2.44 

Baaetnenc 
Vp:  3173  n/s 

Vs:  2723  nVs 

Kp:  a02dB/hi*kHz 
Ks:  a07dBAn*kHz 
p:  173^ 
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APPEHDIZ  D 


A.  MODEL  RESULTS 

Figures  15  through  40  show  the  results  of  the  SAFARI  model 
run  at  a  frequency  of  7  Hz.  OASTL  STLRAN  and  OASTL  SINTGR 
stand  for  Ocean  Acoustic  and  Seismic  Exploration  Synthesis 
Transmission  Loss,  Single  frequency  Transmission  Loss  versus 
Range,  and  Ocean  Acoustic  and  Seismic  Exploration  Synthesis 
Transmission  Loss,  Single  frequency  Integrand.  They  are 
located  on  the  top  right  part  of  the  graphs. 
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o 


(DdL/ZOP)  SS9J^S  PUJJON 


Figura  IS  Transmission  Loss  SAFARI  Float  0 
Cmin«50  m/s,  Cmax>*lE8  m/s,  4096  points,  Ar«7.i41  m. 
Max  FFT  range»29.250  km 


71 


Range  (knn') 


Flgur*  16  Integrand  SAFARI  Float  0 

Cmln«*50  m/s,  Cmax-1E8  m/s,  4096  points,  Ara7.141  m. 
Max  FFT  range*  29.250  km 


rlgur*  17  Transmission  Loss  SAFARI  Float  0 
Cmin«100  m/s,  Cmax»2E8  m/s,  2048  points,  Ar-14.279  m. 
Max  FFT  range«29.243  km 
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rigur*  18  Integrand  SAFARI  Float  0 

Cmln-lOO  m/a,  Cmax«2E8,  2048  points,  Ar»14.279  m, 

Max  FFT  range«29.243  km 
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Flgur*  19  Transmission  Loss  SAFARI  Float  0 
Cmln>740  m/s,  Cmax»lE8  m/s,  4096  points,  Ar«=105.7  m. 
Max  FFT  ranges432.903  km 
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Gi  wavenumber  (10 


Cnin»50  vl/s,  Cmax«lE8  m/s,  4096  points,  Ars7.l41  m, 
Max  FFT  range=29.250  kn 
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Figur*  22  Transmission  Loss  SAFARI  Float  0 
0.25  Contour  Offset 

Cmin»50  m/s,  Cmax»lE8  m/s,  4096  points,  ^>7.141  ro. 
Max  FFT  range«29.250  km 


Figur*  23  Transmission  Loss  SAFARI  Float  0 
No  Contour  Offset 

Cmln«50  m/s,  Cmax-IES  m/s,  4096  points,  Ar«7.141  m. 
Max  FFT  rangeB29.250  km 
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rigur*  24  Transmission  Loss  SAFARI  Float  o 
0.25  Contour  Offset 
Cmin*50  m/s,  Cmax*lE8  m/s,  4096  points,  Ar«7.i41 
Max  FFT  range«29.250  km 


rlgur*  25  Transmission  Loss  SAFARI  Float  0 
0.5  Contour  Offset 

Cmin*50  m/s,  Cmax-1E8  m/s,  4096  points,  ^r»7.l41  m. 
Max  FFT  range«29.250  km 
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Flgur*  27  Transmission  Loss  SAFARI  Float  0 
1.0  Contour  Offset 

Cmin«50  m/s,  Cmax«>lE8  m/s,  4096  points,  iir»7.141  m, 
Max  FFT  range-29.250  km 
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Figure  28  Transmission  Loss  SAFARI  Float  0 
1.25  Contour  Offset 

Cmln«50  m/s,  Cmax«lE8  m/s,  4096  points,  Ars7.l41  m. 
Max  FFT  range»29.250  km 
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riguTtt  29  Transmission  Loss  SAFARI  Float  0 
1.50  Contour  Offset 
Cmin«50  m/s,  cmax«lE8  m/s,  4096  points,  Ar»7.i4l 
Max  FFT  range«29.250  km 


Figur*  30  Transmission  Loss  SAFARI  Float  0 


Standard  Contour  Offset 


Cmin>50  m/s,  Cmax*lE8  m/s,  4096  points,  Ar«7.141  m. 


Max  FFT  range«29.250  km 


* 
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rigtar*  31  Transmission  Loss  SAFARI  Float  0 
Standard  Contour  Offset 

Cmin-50  m/s,  Cmaxa>lE8  m/s,  8192  points,  dr«7.14l  m. 
Max  FFT  range«58.500  km 
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standard  Contour  Offset 


Cinln«740  m/s,  Cmax^ilES  m/B,  4096  points,  Ar=105.7  m, 
Max  FFT  range*432.903  km 
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p6[//QP)  SS9J1S  pLUJON 


Figure  33  Transmission  Loss  SAFARI  Float  0 
Standard  Contour  Offset 

Cmln»740  m/s,  Cmax«lE8  m/s,  4096  points,  Ar-105.7  m. 
Max  FFT  range>432 .903  km 
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Figure  34  Transmission  Loss  SAFARI  Float  0 
Standard  Contour  Offset 

Cmin*740  m/s,  Cmax>=3400  m/s,  4096  points,  Ar=105.7  m. 
Max  FFT  range«432.903  km 
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fe  ♦+0l)  sninpo^^^ 


rlgurs  35  Integrand  SAFARI  Float  0 
Standard  Contour  Offset 

Cmin«740  m/s,  Cmax-3400  m/s,  4096  points,  dr-105. 7  m. 
Max  FFT  range-432.903  km 
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10.0  20.C  30.0  ^0.0  50.0 

Horizontal  wavenumber  (10**  -3) 


standard  Contour  Offset 


Ciiiin«740  m/s,  Cinax-3400  m/a,  1024  points,  dr-105.7  in, 
Max  FFT  range-432.903  km 
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rigiar*  37  integrand  SAFARI  Float  0 
Standard  Contour  Offset 

Cain«740  m/s,  Cmax«3400  m/s,  1024  points,  Ar=105.7  m. 
Max  FFT  range-432.903  km 


Figur*  3F  Transmission  Loss  SAFARI  Float  0 
Standard  Contour  Offset 

Cain«50  m/s,  Cmax-1E8  m/s,  16384  points,  dr«7.141  m. 
Max  FFT  range*!!? . 000  km 


o  o  o  o  p 

o  d  d  d  o 

Th  to  cj 
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Figur*  3f  Integrand  SAFARI  Float  0 
Standard  Contour  Offset 

Cnln*50  n/s,  Cmax«lE8  n/s,  16384  points,  Ar»7.l41  m, 
Max  FFT  range>117.000  km 
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Horizontd  wcvenu.nnber  (10**  —2) 


Figure  40  Transmission  Loss  SAFARI  Float  0 
Water  Column  Only 

Cmin*50  m/s,  Cmax«lE8  m/s,  4096  points,  Ar»7.i41  m. 
Max  FFT  ran9e«29.250  km 
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rigur*  4X  Transmission  Loss,  NATIVE  1,  Event  1,  Float  0, 
7  Hz,  Region  1 
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sstri  uonsfuisuvjx 


Figure  42  Transmission  Loss,  NATIVE  l.  Event  1,  Float  0, 
7  Hz,  Region  2 
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W-> 

f'} 


1  Hz,  Region  1  vs  Region  2 
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Range  (Km) 


I  I  I 


I 


sstri  uoissioisunjx 


Figur*  44  Transmission  Loss,  NATIVE  1,  Event  1,  Float  0, 
10  Hz,  Region  1 
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Ranyc  (ICm) 


SSin  UOISSIUISUBJX 


Flgur*  45  Transmission  Loss,  NATIVE  1,  Event  l,  Float  0, 
10  Hz,  Region  2 
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Figur*  4C  Transmission  Loss,  NATIVE  1,  Event  1,  Float  0, 
10  Hz,  Region  1  vs  Region  2 
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Range  (Km) 


f  S  ^  g  S  S  8 


ss»ri  uoissfiusuRix 


Figur*  47  Transmission  Loss,  SAFARI  vs  NATIVE  1,  Event  1, 
Float  0,  Region  1  at  7  Hz 
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Range  (Km) 


rigur*  48  Transmission  Loss,  SAFARI  vs  NATIVE  1,  Event  1, 
Float  0,  Region  2  at  7  Hz 
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Figurtt  49  Transmission  Loss,  SAFARI  vs  NATIVE  1,  Event  1, 
Float  0,  Region  1  at  7  Hz,  with  2.4  km  Offset  in  Range 
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Range  (Km) 


f  S  ?  R  S  S  S 


swn  UUISStUISUCJX 


rigur*  50  Transmission  Loss,  SAFARI  vs  NATIVE  1,  Event  1, 
Float  0,  Region  2  at  7  Hz,  with  No  Offset  in  Range 
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Range  (Km) 


Flgur*  51  Transmission  Loss,  SAFARI  vs  NATIVE  i,  Event  i, 
Float  0,  Region  1  at  10  Hz 
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Range  (Km) 
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Figur*  S2  Transmission  Loss,  SAFARI  vs  NATIVE  l,  Event  l, 
Float  0,  Region  2  at  10  Hz 
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Range  (Km) 


I 


Figur*  53  Transmission  Loss,  SAFARI  vs  NATIVE  1,  Event  1, 
Float  0,  Region  1  at  10  Hz,  with  0.5  km  Offset  in  Range 
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Figurs  54  Transmission  Loss,  SAFARI  vs  NATIVE  1,  Event  l, 
Float  0,  Region  2  at  10  Hz,  with  -0.75  km  Offset  in  Range 
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